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In  paragraph  2.1,  short-period  and  high-frequency  NORESS  observations 
from  two  aftershock  sequences  offshore  western  Norway,  with  epicentral 
distances  of  424  and  405  km,  respectively,  have  been  used  in  testing  a 
new  hedgehog  inversion  procedure  for  determination  of  source  spectral  scaling 
characteristics,  based  on  a  previously  developed  spectral  ratio  method. 

The  two  main  events  analyzed  have  Hx,  magnitudes  of  about  5  while  the  smallest 
aftershocks  are  about  1.5.  The  observed  spectral  ratios  support  u>2  source 
models  for  both  of  the  earthquake  sequences,  but  only  if  certain  constraints 
are  imposed  on  the  scaling  parameter  that  determines  the  stress  drop  sensitivity 
with  respect  to  seismic  moment. 

P-wave  focusing  effects  at  the  NORSAR  array  for  Novaya  Zemlya  explosions 
are  investigated  in  paragraph  2.2.  It  is  found  that  the  amplitude  pattern 
for  Novaya  Zemlya  explosions  shows  an  order  of  magnitude  variation  across 
the  60  km  aperture  NORSAR  array.  The  seismometer  at  site  03C01  records 
the  strongest  signals,  and  it  is  Bhown  by  simple  SNR  scaling  that  this 
seismometer  would  record  detectable  signals  for  Novaya  Zemlya  explosions 
of  nn  =  2.0 — 2.5  during  normal  noise  conditions. 

A  major  consideration  in  selecting  the  NORESS  array  site  has  been  the 
high  amplification  of  short-period  P  signals  from  directions  to  the  east. 

An  investigation,  described  in  paragraph  2.3,  has  been  initiated  aimed 
at  characterizing  the  structure  responsible  for  the  presumed  wave  focusing. 

The  data  base  provided  by  NORESS  itself  must  be  extended  for  this  purpose, 
and  we  have  started  processing  data  from  the  large-aperture  NORSAR  array 
that  encompasses  the  NORESS  site.  At  this  stage  a  method  has  been  developed 
to  downward  extrapolate  the  observed  surface  wavefield  at  NORSAR.  The  method 
is  posed  as  an  inverse  problem  taking  into  account  the  limited  areal  extent 
of  sampling  at  the  surface.  Preliminary  results  suggest  that  the  source 
of  focusing  at  NORESS  is  far  below  the  Moho. 

An  important  incentive  for  investigating  the  natural  seismicity  of  Norway 
and  surrounding  areas  is  that  the  nuclear  explosion  monitoring  programs 
need  the  best  possible  background  information  about  the  seismicity,  not 
only  in  the  areas  where  such  explosions  could  take  place,  but  also  in  the 
areas  within  local  and  regional  distance  from  the  monitoring  stations. 

In  the  new  Intelligent  Honitoring  System  now  installed  at  NORSAR,  such 
background  seismicity  can  be  displayed  as  overlays.  For  this  purpose, 
three  different  catalogues  have  been  prepared,  covering  the  three  time 
periods  1880 — 1954,  1955 — 1979  and  1980 — 1989.  The  purpose  of  the  contribution 
in  paragraph  2.4  is  to  document  these  catalogues,  and  to  outline  briefly 
their  seismotectonic  context. 

Paragraph  2.5  contains  a  report  on  the  scientific  symposium  entitled 
*  ‘Regional  Seismic  Arrays  and  Nuclear  Test  Ban  Verification”,  held  in 
Oslo,  Norway,  during  14 — 17  February  1990.  The  purpose  of  the  symposium 


was  to  assess  the  state-of-the-art  research  on  regional  seismic  arrays 

and  associated  topics.  In  particular,  the  symposium  focused  upon  the  advanced 

regional  arrays  NORESS  and  ARCESS  in  Norway  and  their  associated  data  processing 

facilities,  in  the  light  of  the  potential  of  such  arrays  to  provide  a  much 

improved  monitoring  capability  for  a  future  comprehensive  nuclear  test 

ban  treaty.  In  paragraph  2.5,  we  give  a  brief  review  of  some  of  the  results 

presented  during  the  symposium.  The  majority  of  the  papers  presented  will 

be  published  in  a  special  issue  of  the  Bulletin  of  the  Seismological  Society 

of  America,  scheduled  to  appear  in  December  1990.  Authors  and  titles  of 

the  28  papers  of  the  special  issue  are  given. 
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Under  Contract  No.  F49620-C-89-0038,  NTNF/NORSAR  is  conducting 
research  within  a  wide  range  of  subjects  relevant  to  seismic  monitoring.  The 
emphasis  of  the  research  program  is  on  developing  and  assessing  methods 
for  processing  of  data  recorded  by  networks  of  small-aperture  arrays  and  3- 
component  stations,  for  events  both  at  regional  and  teleseismic  distances.  In 
addition,  more  general  seismological  research  topics  are  discussed. 

Each  quarterly  technical  report  under  this  contract  presents  one  or  several 
separate  investigations  addressing  specific  problems  within  the  scope  of  the 
statement  of  work.  Summaries  of  the  research  efforts  within  the  program  as  a 
whole  are  given  in  annual  technical  reports. 

This  Scientific  Report  No.  7  is  the  annual  technical  report  for  the  period  1 
October  1989  -  30  September  1990.  It  contains  five  separate  contributions,  and 
also  abstracts  for  the  investigations  submitted  as  quarterly  technical  reports 
during  FY90. 
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1  Summary 


This  Annual  Technical  Report  describes  the  work  accomplished  under  Con¬ 
tract  No.  F49620-C-89-0038  during  the  period  1  October  1989  -  30  September 
1990.  The  report  contains  five  separate  contributions  (paragraphs  2.1  through 
2.5),  and  in  addition  abstracts  of  the  investigations  submitted  as  quarterly 
technical  reports  during  FY90  (paragraph  2.6). 

In  paragraph  2.1,  short-period  and  high-frequency  NORESS  observations 
from  two  aftershock  sequences  offshore  western  Norway,  with  epicentral  dis¬ 
tances  of  424  and  405  km,  respectively,  have  been  used  in  testing  a  new  hedge¬ 
hog  inversion  procedure  for  determination  of  source  spectral  scaling  character¬ 
istics,  based  on  a  previously  developed  spectral  ratio  method.  The  two  main 
events  analyzed  have  Mz  magnitudes  of  about  5  while  the  smallest  aftershocks 
are  about  1.5.  The  observed  spectral  ratios  support  u2  source  models  for  both 
of  the  earthquake  sequences,  but  only  if  certain  constraints  are  imposed  on 
the  scaling  parameter  that  determines  the  stress  drop  sensitivity  with  respect 
to  seismic  moment. 

P-wave  focusing  effects  at  the  NORSAR  array  for  Novaya  Zemlya  explo¬ 
sions  are  investigated  in  paragraph  2.2.  It  is  found  that  the  amplitude  pattern 
for  Novaya  Zemlya  explosions  shows  an  order  of  magnitude  variation  across 
the  60  km  aperture  NORSAR  array.  The  seismometer  at  site  03C01  records 
the  strongest  signals,  and  it  is  shown  by  simple  SNR  scaling  that  this  seis¬ 
mometer  would  record  detectable  signals  for  Novaya  Zemlya  explosions  of  mt 
=  2.0-2.5  during  normal  noise  conditions. 

A  major  consideration  in  selecting  the  NORESS  array  site  has  been  the 
high  amplification  of  short-period  P  signals  from  directions  to  the  east.  An 
investigation,  described  in  paragraph  2.3,  has  been  initiated  aimed  at  charac¬ 
terizing  the  structure  responsible  for  the  presumed  wave  focusing.  The  data 
base  provided  by  NORESS  itself  must  be  extended  for  this  purpose,  and  we 
have  started  processing  data  from  the  large-aperture  NORSAR  array  that 
encompasses  the  NORESS  site.  At  this  stage  a  method  has  been  developed 
to  downward  extrapolate  the  observed  surface  wavefield  at  NORSAR.  The 
method  is  posed  as  an  inverse  problem  taking  into  account  the  limited  areal 
extent  of  sampling  at  the  surface.  Preliminary  results  suggest  that  the  source 
of  focusing  at  NORESS  is  far  below  the  Moho. 

An  important  incentive  for  investigating  the  natural  seismicity  of  Norway 
and  surrounding  areas  is  that  the  nuclear  explosion  monitoring  programs  need 
the  best  possible  background  information  about  the  seismicity,  not  only  in  the 
areas  where  such  explosions  could  take  place,  but  also  in  the  areas  within  local 
and  regional  distance  from  the  monitoring  stations.  In  the  new  Intelligent 
Monitoring  System  now  installed  at  NORSAR,  such  background  seismicity 
can  be  displayed  as  overlays.  For  this  purpose,  three  different  catalogues  have 
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been  prepared,  covering  the  three  time  periods  1880-1954,  1955-1979  and 
1980-1989.  The  purpose  of  the  contribution  in  paragraph  2.4  is  to  document 
these  catalogues,  and  to  outline  briefly  their  seismotectonic  context. 

Paragraph  2.5  contains  a  report  on  the  scientific  symposium  entitled  “Re¬ 
gional  Seismic  Arrays  and  Nuclear  Test  Ban  Verification”,  held  in  Oslo,  Nor¬ 
way,  during  14-17  February  1990.  The  purpose  of  the  symposium  was  to 
assess  the  state-of-the-art  research  on  regional  seismic  arrays  and  associated 
topics.  In  particular,  the  symposium  focused  upon  the  advanced  regional  ar¬ 
rays  NORESS  and  ARCESS  in  Norway  and  their  associated  data  processing 
facilities,  in  the  light  of  the  potential  of  such  arrays  to  provide  a  much  improved 
monitoring  capability  for  a  future  comprehensive  nuclear  test  ban  treaty.  In 
paragraph  2.5,  we  give  a  brief  review  of  some  of  the  results  presented  during 
the  symposium.  The  majority  of  the  papers  presented  will  be  published  in  a 
special  issue  of  the  Bulletin  of  the  Seismological  Society  of  America,  scheduled 
to  appear  in  December  1990.  Authors  and  titles  of  the  28  papers  of  the  special 
issue  are  given. 
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2  Summary  of  Technical  Findings  and  Accomplish¬ 
ments 

2.1  Source  spectral  scaling  inversion  for  two  earthquake  se¬ 
quences  offshore  western  Norway 

Introduction 

The  way  in  which  earthquake  source  spectral  characteristics  depend  on 
seismic  moment  (or  magnitude)  is  important  in  the  many  situations  in  which 
observations  from  earthquakes  in  a  certain  magnitude  range  are  used  for  pre¬ 
dicting  effects  from  earthquakes  in  another  (usually  higher)  magnitude  range. 
The  discussion  of  source  scaling  laws  often  focuses  on  the  self-similarity  prin¬ 
ciple  (Aki,  1967),  a  departure  from  which  is  frequently  seen  in  terms  of  non¬ 
constant  stress  drop  (with  seismic  moment),  and  sometimes  also  in  terms  of 
non-exponential  frequency-magnitude  distributions  (Aki,  1987), 

The  departure  from  the  constant  stress  drop  (CSD)  situation  invariably 
means  increasing  stress  drop  (ISD)  with  seismic  moment,  corresponding  to  a 
reduced  rate  of  increase  of  corner  frequencies  and  a  reduced  rate  of  decrease 
of  source  radii  as  magnitude  decreases.  While  the  large  number  of  reports  on 
such  observations  often  agree  that  they  start  to  occur  at  a  magnitude  of  around 

3  and/or  a  moment  of  around  1021  Nm,  similar  agreements  are  less  common 
with  regard  to  explanation  (Dyshart  et  al .,  1988).  These  fall  generally  in  two 
categories,  propagation  and/or  site  effects,  and  source  effects. 

A  central  concept  in  the  first  case  is  the  fmax  (Hanks,  1982),  often  but  not 
always  interpreted  as  a  site  effect  parameter,  appearing  for  large  earthquakes 
as  a  second  corner  frequency  and  for  smaller  events  as  a  limiting  value  for 
the  corner  frequency.  The  second  type  of  explanations  may  involve  concepts 
such  as  characteristic  barrier  spacings  and  critical  minimum  crack  length,  but 
the  reason  for  the  departure  from  self-similarity  could  also  simply  be  that  the 
corner  frequency  of  small  earthquakes  are  controlled  by  the  fault  zone  width 
rather  than  by  the  fault  length  (Aki,  1987). 

Spatially  confined  earthquake  sequences  have  been  used  successfully  in  de¬ 
termination  of  source  spectral  characteristics.  While  earlier  studies  (Chael, 
1987;  Chael  and  Kromer;  1988;  Chun  et  al .,  1989)  only  compared  spectral 
ratio  curve  shapes  to  theoretical  ones,  Boore  and  Atkinson  (1989)  extended 
the  method  by  using  independently  determined  seismic  moments,  thereby  con¬ 
trolling  also  the  absolute  levels  of  the  spectral  ratios.  These  studies  all  find 
support  for  an  u/2  high-frequency  fall-off,  while  stress  drop  (when  determined) 
follows  both  CSD  and  ISD  models. 

The  data  used  in  the  present  paper  coincide  in  part  with  those  used  by 
Chael  and  Kromer  (1988),  who  analyzed  an  1986  aftershock  sequence  off  the 
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western  coast  of  Norway.  In  addition,  we  use  data  also  from  an  1989  after¬ 
shock  sequence  not  far  from  the  one  in  1986.  We  use  essentially  the  same 
spectral  ratio  technique,  but  depart  from  earlier  studies  in  that  we  replace 
their  trial-and-error  method  with  one  in  which  a  hedgehog  inversion  technique 
is  applied.  This  provides  (given  that  the  data  are  of  sufficient  quality)  more  re¬ 
liable  solutions,  and  better  possibilities  for  analyzing  resolution,  uncertainties, 
and  trade-offs  between  parameters. 

Data 

The  data  used  in  this  analysis  comprise  high-frequency  (hf)  and  short- 
period  (sp)  recordings  at  one  of  the  sites  in  the  NORESS  regional  array  in 
Norway  (Mykkeltveit  et  a/.,  1990),  from  two  aftershock  sequences  offshore 
western  Norway  in  1986  and  1989.  The  two  sequences,  here  named  A  and 
B ,  are  identified  in  Table  2.1.1,  with  origin  times,  magnitudes  and  seismic 
moments.  Aq  and  Bo  refers  to  the  main  events,  and  we  have  analyzed  six 
aftershocks  in  each  case,  At  -  Aq  and  2?i  -  Be-  It  is  seen  from  Table  2.1.1  that 
the  magnitudes  of  the  main  events  are  4.8  and  4.9,  while  the  aftershocks  are 
all  in  the  range  2.8  to  1.4.  For  sequence  B  all  of  the  aftershocks  occur  within 
3  days,  while  sequence  A  is  scattered  over  a  time  period  of  8  days. 

Event  Aq  occurred  off  the  west  coast  of  Norway  at  62.71°N,  4.69°E,  at 
a  distance  of  424  km  from  the  NORESS  array,  and  Bo  occurred  at  61.97°N, 
4.42°E  at  a  distance  of  405  km.  The  distance  between  the  event  sequences  is  84 
km,  both  of  them  occurred  on  the  landward  side  of  the  continental  margin  (B 
closest  to  the  coast),  and  the  focal  mechanisms  in  both  cases  indicate  thrust 
faulting  in  response  to  EW  compressional  stress.  More  detailed  descriptions 
of  the  events  are  given  by  Hansen  et  al.  (1989),  and  for  more  details  about 
the  seismotectonics  of  the  region  we  refer  to  Bungum  et  al.  (1991).  There  are 
no  indications  of  significantly  different  tectonic  conditions  between  the  two 
epicentral  areas.  Event  3o  has  been  located,  using  near-source  recordings, 
very  close  to  (possibly  below)  the  crust-mantle  boundary  (24-28  km),  while  A0 
most  probably  is  somewhat  shallower  (Hansen  et  al.,  1989).  In  the  subsequent 
analysis,  only  vertical  component  data  have  been  used. 

For  the  aftershock  data  we  have,  in  addition  to  the  high-frequency  (hf) 
data,  also  used  short-period  (sp)  data  from  instruments  located  very  close  to 
the  hf  elements,  partly  for  the  purpose  of  stabilizing  the  estimates  and  partly 
in  order  to  provide  records  also  for  aftershocks  As  and  A$,  for  which  hf  records 
where  missing  (see  Table  2.2.1).  The  dynamically  more  limited  sp  data  are 
reaching  clipping  level  for  Lg,  but  could  still  be  used  because  we  analyze  data 
only  between  the  Pn  and  the  5n  arrivals. 

Method  of  Analysis 

The  main  advantage  of  spectral  ratio  methods,  which  are  used  in  many 
different  types  of  applications,  is  that  both  site  effects  (including  instrument 
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response)  and  path  effects  (earth  transfer  functions)  are  removed  in  taking 
the  ratios.  Under  certain  conditions,  source  radiation  effects  are  also  removed, 
while  earth  noise  effects  are  additive  and  should  be  dealt  with  specifically. 

The  spectral  ratio  method  has  been  used  recently  in  source  scaling  inves¬ 
tigations  by  Chael  (1987),  Chael  and  Kromer  (1988),  Chun  et  al.  (1989)  and 
Boore  and  Atkinson  (1989),  and  we  refer  to  these  papers  for  a  more  detailed 
description.  The  relationship  for  the  power  spectral  ratio  R2(d)  between  the 
i’th  event  (aftershock)  and  the  reference  event  (main  shock)  r  in  this  case  is 


where  u  is  angular  frequency,  u>q  is  corner  frequency,  and  Mq  is  seismic  mo¬ 
ment.  The  spectral  ratios  are  seen  to  be  controlled  here  by  the  two  parameters 
7  and  6,  where  7  controls  the  shape  of  the  source  spectrum  as  demonstrated 
by  the  following  relationship  (Aki,  1967;  Aki  and  Richards,  1980)  between  the 
far-field  displacement  spectrum  Q(w)  and  its  low-frequency  asymptote  fl(0) 
(which  in  turn  scales  linearly  with  seismic  moment) 


Values  of  2  and  3  for  7  correspond  to  the  w2  and  the  w3  models,  respectively. 
The  parameter  £  in  equation  (1)  controls  the  scaling  relation  between  corner 
frequencies  and  seismic  moment  as  follows 

MqUqS  =  constant  (3) 


/1  _  3  \ 

which  in  turn  leads  to  a  proportionality  between  stress  drop  Act  and  Mq  , 
or  alternatively 

=  A  (4) 


where  subscripts  i  and  r  are  defined  as  for  equation  (1).  Values  of  3  and  4  for 
6  correspond  to  CSD  and  ISD  (Nuttli,  1983a)  models,  respectively.  Unless  an 
independently  determined  estimate  of  the  stress  drop  for  the  reference  event 
is  available,  the  spectral  ratio  method  only  provides  the  stress  drop  ratios. 


A  commonly  used  realization  of  the  suite  of  source  models  allowed  by  the 
different  combinations  of  7  and  8  is  the  one  due  to  Brune  (1970;  1971),  which 
is  an  w2  CSD  model.  There  are  many  other  models  that  also  relate  corner 
frequencies  to  stress  drop. 
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Data  Analysis 

The  data  have  been  analyzed  in  this  study  essentially  as  described  in  de¬ 
tail  by  Chael  (1987),  Chael  and  Kromer  (1988)  and  Chun  et  al.  (1989).  We 
have  used  the  data  between  the  arrivad  of  the  Pn  phise  and  the  arrival  of 
the  Sn  phase,  and  computed  power  spectra  using  the  direct  method  of  power 
spectral  estimation  (the  Fourier  transform  of  the  data  multiplied  by  its  com¬ 
plex  conjugate),  combined  with  a  block- averaging  with  overlapping  blocks  for 
the  purpose  of  spectral  smoothing.  The  approximate  stationarity  of  the  P 
coda  makes  it  possible  to  estimate  power  spectra  in  this  way,  and  the  block¬ 
averaging  is  simply  one  of  the  ways  in  which  spectral  stability  is  achieved  at 
the  expense  of  frequency  resolution. 

For  both  of  the  events,  about  35  sec  of  data  were  available  between  the 
Pn  and  the  Sn  phases.  Signal  spectra  have  been  estimated,  using  the  block¬ 
averaging  method,  for  each  of  the  10  hf  and  12  sp  records  (see  Table  2.1.1), 
and  additional  smoothing  has  been  achieved  in  the  frequency  domain.  The 
rationale  for  smoothing  the  spectra  so  extensively  is  that  dynamic  stability  in 
the  present  analysis  is  more  important  than  frequency  resolution. 

Noise  spectra  have  been  estimated  in  the  same  way,  and  for  equally  long 
time  windows  preceding  the  Pn  arrivals  (except  for  event  B\  where  only  10 
sec  of  noise  were  available).  These  noise  spectra  are  necessary  for  two  rea¬ 
sons,  because  noise  power  should  be  subtracted  from  signal  power  (se  detailed 
discussion  on  this  subject  by  Chun  et  al .,  1989)  and  because  a  noise  estimate 
is  needed  in  order  to  decide  which  parts  of  the  noise-corrected  spectrum  that 
can  be  used  in  each  particular  case. 

The  signal- to- noise  ratios  (SNR’s)  for  the  main  events  were  around  30  dB 
for  .4o  and  70  dB  for  Bo,  while  the  aftershocks  had  SNR  values  of  up  to  11 
dB  for  sequence  A  and  17  dB  for  B.  There  were  also,  however,  some  low- 
frequency  SNR  values  in  both  cases  at  or  below  zero,  a  finding  which  has  been 
used  in  determining  the  exact  frequency  band  of  observations  to  use  in  each 
particular  case.  We  have  found  that  subtracting  the  noise  makes  a  difference 
only  for  the  weakest  ones  of  the  aftershocks. 

In  earlier  analyses  of  source  scaling  characteristics  using  the  spectral  ratio 
method  it  has  been  common  to  compare  observations  only  to  a  suite  of  theo¬ 
retical  spectral  ratio  curves,  for  the  purpose  of  comparing  spectral  shapes.  It 
is  obvious,  however,  that  the  availability  of  independently  determined  seismic 
moments  can  assist  in  constraining  the  suite  of  acceptable  models,  as  pointed 
out  also  by  Boore  and  Atkinson  (1989).  We  therefore  proceed  now  with  estab¬ 
lishing  seismic  moments  for  the  data  analyzed  in  this  study,  before  starting 
the  source  scaling  inversion  analysis. 

The  magnitudes  used  in  the  present  paper  have  been  estimated  using  a 
magnitude  inversion  algorithm  recently  applied  by  Alsaker  et  al.  (1991),  who 


used  a  data  base  of  741  records  from  195  earthquakes  in  developing  attenua¬ 
tion  terms  and  station  corrections.  The  resulting  Mi  magnitudes  are  on  the 
average  about  0.4  units  below  those  obtained  with  the  earlier  scale  (Bath  et 
al.,  1976).  When  applied  to  the  data  analyzed  here,  the  new  scale  gives  Mi 
values  as  shown  in  Table  2.1.1.  While  only  NOR.ESS  data  were  available  for 
magnitude  assessment  for  sequence  A,  we  could  for  the  B  events  use  up  to  6 
stations  as  shown  in  Table  2.1.1. 

Seismic  moments  were  estimated  following  an  inversion  algorithm  recently 
applied  by  Kvamme  and  Hansen  (1991)  to  essentially  the  same  data  base  that 
Alsaker  et  al.  (1991)  used  for  the  magnitude  scale.  The  inversion  algorithm 
(see  also  Sereno  et  al .,  1988)  utilizes  ground  motion  displacement  amplitude 
spectra,  and  we  have  tested  it  here  with  different  options  (for  which  parameters 
to  determine,  always  including  moments)  and  different  data  bases  (always 
including  the  two  aftershock  sequences)  in  order  to  assess  the  stability  of 
the  results.  The  seismic  moments  determined  are  also  in  this  case  shown 
in  Table  2.1.1,  and  also  this  time  there  are  many  more  records  available  for 
the  B  events.  The  additional  station  for  event  Ao  in  this  case  is  a  NOR.ESS 
broadband  (intermediate  period,  ip)  record  that  helps  to  stabilize  the  moment 
estimates  for  the  main  events  The  aftershocks  were  too  weak  to  be  recorded 
on  the  ip  sensors. 

The  spectral  data  for  the  two  earthquake  sequences  have,  together  with 
the  seismic  moments,  been  used  in  a  simple  hedgehog  inversion  procedure  for 
the  purpose  of  determining  simultaneously  the  scaling  parameters  7,  6  and 
wo,r>  with  and  without  a  redetermination  of  the  moment  ratios.  The  results 
of  this  analysis  are  shown  in  Table  2.1.2  and  in  Figures  2.1.1  and  2.1.2,  for 
sequences  A  and  B,  respectively.  It  is  seen  there  that  6  is  3.24  for  sequence 
A  (close  to  an  CSD  model)  while  it  is  3.80  for  sequence  B  (close  to  an  ISD 
model).  We  will  return  to  these  results  after  discussing  the  method  in  more 
detail. 

The  starting  point  for  the  inversion  is  equation  (1),  which  allows  us  to 
estimate  the  spectral  ratios  as  functions  of  frequency,  given  7,  6 ,  u0yT  and 
the  moment  ratios.  The  problems  with  linearizing  equation  (1)  have  been 
circumvented  here  by  using  a  hedgehog  approach  (Aki  and  Richards,  1980), 
which  in  fact  is  particularly  appropriate  in  this  case  because  we  have  physical 
limitations  that  helps  us  to  confine  the  parameter  space  to  within  certain 
limits.  For  the  parameter  7,  it  suffices  initially  to  consider  only  the  values 
of  2  (w2-model)  and  3  (w3-model),  the  parameter  6  can  be  limited  to  values 
centered  in  the  range  3  (CSD)  to  4  (ISD),  and  the  corner  frequency  /o,r  (= 
u>o,r/2jr)  of  the  main  event  is  bound  by  the  estimated  seismic  moment.  Seismic 
moments  of  the  order  of  1016  Nm  should  correspond  to  corner  frequencies  of 
the  order  of  1-2  Hz. 

Given  these  parameter  space  constraints,  we  computed  initially  theoretical 
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spectral  ratios  using  equation  (1)  for  a  sufficiently  dense  grid  of  parameter 
combinations  within  the  limits  2-3  for  7, 2.5-4.5  for  6,  and  0.5-3.0  for  /o,r.  For 
each  parameter  combination  we  computed  the  root  mean  square  (RMS)  of  the 
deviations  between  observed  and  computed  log-amplitudes,  averaged  over  all 
frequencies  (with  acceptable  spectral  SNR)  and  over  all  of  the  aftershocks.  A 
proper  weighting  in  frequency  domain  was  obtained  by  resampling  the  spectra 
to  an  equidistant  spacing  in  log-frequency  (20  samples  per  decade)  before  the 
RMS  computation  (by  interpolating  at  low  frequencies  and  smoothing  at  high 
frequencies),  while  a  similar  weighting  between  the  events  was  obtained  by 
using  log-amplitudes.  In  the  cases  when  both  hf  and  sp  data  were  available  a 
weight  of  0.5  was  given  to  each  of  them.  In  this  way,  the  resulting  (properly 
normalized)  RMS  values,  as  computed  from  log-amplitude  differences,  can  be 
related  directly  and  easily  to  average  amplitude  ratios. 

This  search  in  parameter  space  can  either  be  conducted  with  a  sufficiently 
dense  grid  initially,  or  one  can  use  a  coarser  grid,  determine  a  minimum,  and 
iterate  with  gradually  finer  grids.  The  difference  between  these  two  approaches 
is  only  computer  time.  We  have  also,  however,  developed  an  alternative  option 
under  which  we  for  each  grid  point  also  search  for  the  moment  ratio  that  gives 
the  best  fit  between  computed  and  observed  spectral  ratios.  The  moment  ratio 
corresponding  to  the  best-fitting  parameter  combination  in  the  first  iteration 
is  then  carried  into  the  next  iteration,  where  the  grid  is  a  bit  narrower,  and 
the  procedure  continues  until  a  certain  convergence  criterion  is  fulfilled  (we 
have  used  1%  change).  In  that  case  the  iterations  are  essential  for  the  results 
because  of  the  gradual  adjustments  of  the  moment  ratios.  It  is  this  latter 
option  that  has  been  used  when  obtaining  the  results  shown  in  Figures  2.1.1 
and  2.1.2. 

When  testing  the  inversion  procedure  it  soon  became  clear,  when  confining 
the  parameters  to  the  initially  defined  range,  that  an  w2-model  satisfied  the 
observations  much  better  than  an  w3- model,  as  found  also  by  most  others  using 
the  spectral  ratio  method.  We  will  return  to  the  question  of  the  7  value  later, 
while  in  the  following  presenting  the  results  for  a  fixed  7  of  2.  Figure  2.1.3 
shows,  for  aftershock  sequence  B,  the  results  for  the  first  four  iterations  when 
inverting  also  for  moment.  The  two  parameters  to  be  resolved  are  6  and  /o,r, 
and  contour  diagrams  are  used  to  depict  the  RMS  levels.  There  is  one  more 
iteration  in  each  of  these  cases,  but  without  any  real  change  in  the  location  of 
the  minima. 

The  first  impression  we  get  from  Figure  2.1.3  is  the  elongation  of  the  con¬ 
tours,  indicating  that  the  resolution  is  non-uniform.  We  notice  moreover  only 
a  very  small  shift  in  the  locations  of  the  minima  over  the  iterations,  with  <5//o,r 
changing  from  3.70/1.50  to  3.80/1.53.  The  results  are  similarly  stable  for  se¬ 
quence  A.  From  Figure  2.1.3  we  also  find  that  the  inversion  for  moment  ratios 
does  not  really  influence  the  results  for  the  other  parameters  determined  in 
the  inversion. 
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The  results  from  the  present  source  parameter  inversion  are  summarized 
in  Table  2.1.2,  where  the  new  moment  values  already  have  been  discussed.  For 
the  aftershock  corner  frequencies,  equation  (3)  gives  us  the  relation 


while  equation  (4)  has  been  used  in  computing  the  stress  drop  ratios.  An 
example  of  what  this  could  give  in  terms  of  stress  drops  is  shown  in  the  last 
column  of  Table  2.1.2,  where  a  Brune  model  (equation  (5))  has  been  used  in 
computing  reference  event  stress  drop.  It  is  seen  clearly  there  that  a  8- value  of 
3.24  for  sequence  A  gives  only  moderately  varying  stress  drops  with  moment 
(from  100  to  53%),  while  a  6- value  of  3.80  for  sequence  B  gives  a  much  stronger 
stress  drop  variation  (from  100  to  10%).  As  noted  earlier,  a  CSD  (constant 
stress  drop)  model  corresponds  to  6=  3.0  and  an  ISD  (increasing  stress  drop) 
model  to  6=  4.0. 

It  was  mentioned  above  that  a  number  of  options  were  considered  in  the 
initial  determination  of  seismic  moments  (Kvamme  and  Hansen,  1991),  and  we 
extended  this  also  to  testing  a  number  of  the  resulting  solutions  in  the  present 
inversion.  The  results  were  found  to  be  quite  stable  for  all  of  the  parameters 
except  the  main  event  corner  frequency,  where  significant  changes  were  found 
by  using  different  (but  still  reasonable)  sets  of  initial  moments.  Naturally,  the 
seismic  moment  for  the  main  event  was  most  critical  in  this  respect.  While 
this  means  that  the  corner  frequency  resolution  is  relatively  poor,  it  is  still 
important  to  include  this  parameter  in  the  inversion  since  the  use  of  a  fixed 
value  otherwise  easily  could  bias  the  results  for  the  other  parameters. 

It  was  noted  above  that  the  initial  solutions  were  much  better  for  an  w2 
source  model  as  compared  to  an  w3-model.  In  noting,  however,  that  the 
best  fit  in  the  latter  case  was  obtained  for  simultaneously  high  values  of  6 
and  /o,r,  we  extended  the  ranges  of  these  parameters  beyond  their  initial 
limits.  What  we  found  then  were  the  results  shown  in  Figure  2.1.4  (sequence 
B  only),  where  the  box  in  the  lower-left  corner  indicates  the  previously  used 
parameter  range.  The  results  are  reasonably  similar  also  for  sequence  A,  and 
we  find  again  that  the  moment-inversion  option  does  not  change  the  results 
significantly  for  either  of  the  aftershock  sequences.  The  fit  between  computed 
and  observed  (in  terms  of  RMS)  is  almost  as  good  as  earlier,  which  means  that 
there  is  a  basic  ambiguity  between  the  u>2  and  the  u3  source  models  that  can 
be  resolved  only  by  invoking  independent  physical  constraints  for  S  and/or  /o, 
or  by  determining  7  independently  from  near-source  observations. 

When  considering  only  the  initial  parameter  range  in  Figure  2.1.4  we  see 
that  the  RMS  values  are  mostly  above  1.0,  as  compared  to  values  around  0.20 
when  using  an  w2  source  model.  This  was  the  basis  for  initially  rejecting  the 
u3  model.  Comparisons  between  computed  and  observed  spectral  ratios  (as 
in  Figures  2.1.1  and  2.1.2)  also  indicate  a  significantly  poorer  fit  in  this  case. 
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The  good  fit  for  the  w3-based  solutions  in  Figure  2.1.4  makes  it  desirable 
to  investigate  this  finding  further,  because  of  its  potential  importance  with 
respect  to  trade-off  effects  between  parameters.  We  did  this  by  looping  the 
source  spectral  slope  parameter  7  between  values  of  1.5  and  3.0,  with  results 
as  given  in  Figure  2.1.5,  for  both  of  the  earthquake  sequences.  What  we  see 
there  is  that  each  of  these  7  values  in  fact  has  one  combination  of  6  and  /o 
that  gives  a  very  good  fit  between  computed  and  observed  spectral  ratios, 
with  corresponding  RMS  values  in  the  range  0.18-0.21  for  sequence  A  and 
0.14-0.21  for  B.  Even  though  there  for  sequence  A  is  a  weak  RMS  minimum 
for  the  more  commonly  accepted  6  values  (between  the  dotted  lines)  we  do  not 
attribute  any  real  significance  to  this,  especially  when  considering  that  this 
sequence  has  less  resolving  power  than  B.  It  is  noteworthy,  however,  that  the 
scaling  parameter  6  varies  quite  smoothly  and  systematically  with  7,  which  in 
turn  increases  the  confidence  of  the  estimated  6  values  of  3.24  and  3.80  for  the 
two  sequences,  given  a  7  of  2.0. 

Discussion 

The  results  obtained  in  this  study  are  in  general  consistent  with  results 
from  earlier  studies  using  the  same  method.  Chael  and  Kromer  (1988),  for 
example,  analyzing  also  earthquake  sequence  A,  found  that  an  u;3  source  model 
and  a  7  of  3  (CSD)  fitted  their  observations  best.  For  an  w3  source  model  we 
have  found  a  7  of  3.24  for  the  same  earthquake  sequence.  The  use  of  seismic 
moments  in  this  study  makes  it  more  interesting  to  compare  results  with  those 
of  Boore  and  Atkinson  (1989),  however,  who  emphasized  the  non-uniqueness 
in  the  solutions  if  only  spectral  ratios  are  considered.  They  found,  for  example, 
that  good  spectral  ratio  fits  could  be  maintained  by  simultaneously  increasing 
6  and  decreasing  the  stress  drop  (/0  in  the  present  case),  a  finding  which 
corresponds  to  following  the  elongated  ridges  in  Figure  2.1.3.  Also,  when 
comparing  w3  to  w3  models  they  found  that  a  much  higher  6  was  needed  in 
the  latter  case,  which  corresponds  to  what  can  be  seen  from  Figure  2.1.4. 

Those  observations  can  serve  as  a  useful  prelude  to  one  of  the  main  results 
of  this  study,  namely  the  basic  non-uniqueness  demonstrated  in  Figure  2.1.5. 
There,  a  very  large  range  of  parameter  values  is  used,  and  still  found  to  be 
consistent  with  observed  spectral  ratios.  In  that  case  it  looks  as  if  we  have 
practically  no  resolving  power  with  respect  to  the  scaling  parameter  7. 

We  can,  as  mentioned  earlier,  choose  between  these  solutions  only  by  invok¬ 
ing  external  constraints  on  the  derived  parameters.  This  applies  in  particular 
to  6,  where  a  restriction  to  the  commonly  accepted  range  between  3  (CSD) 
and  4  (ISD)  immediately  restricts  7  to  within  the  range  1. 9-2.3  for  earthquake 
sequence  A  and  1.7-2. 1  for  B  (see  Figure  2.1.5).  It  is  important  to  note  in  this 
respect,  however,  that  McGarr  (1986)  has  found  6  values  of  up  to  above  6  for 
earthquake  sequences  in  the  same  magnitude  range  as  ours,  with  indications 
that  the  largest  ones  within  at  least  for  some  of  the  sequences  scale  according 
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to  a  CSD  model  of  3).  The  transition  in  6  was  found  to  be  in  the  range  1014 
to  1017  Nm,  depending  on  tectonic  regime  (McGarr,  1984).  It  is  therefore  not 
possible  to  exclude  per  se  the  higher  6  values  in  Figure  2.1.5. 

Given  this  non-unique  situation,  the  range  of  acceptable  solutions  could 
of  course  also  be  limited  if  sufficiently  reliable  near-source  observations  had 
been  available  to  allow  for  a  direct  estimation  of  7  (which  more  often  than  not 
results  in  u>2  models).  Unfortunately,  no  such  data  are  available  in  the  present 
case.  The  corner  frequency  also  enters  this  discussion,  however,  since  we  see 
from  Figure  2.1.4  that  7  =  3  gives  about  50%  higher  /o  values  (3.20  and  2.25 
Hz)  as  7  =  2  (2.19  and  1.53  Hz),  which  in  turn  are  somewhat  higher  than 
those  obtained  from  the  inversion  algorithm  of  Kvamme  and  Hansen  (1991). 
We  therefore  have  some,  albeit  not  strong,  evidence  for  choosing  the  solutions 
resulting  from  7  =  2  (and  6  =  3-4).  The  use  of  this  constraint  on  7  should  be 
kept  in  mind  in  the  following. 

Earthquake  sequence  B  has  better  quality  and  more  resolving  power  than 
A,  which  covers  a  smaller  magnitude  range.  An  important  additional  point  in 
favor  of  B  is  that  many  more  seismic  records  than  those  used  in  computing 
spectral  ratios  in  that  case  have  been  used  in  computing  magnitude  and  seismic 
moments.  Another  point  related  to  the  question  of  data  quality  is  that  the 
spectral  ratio  method  ideally  requires  the  same  epicenter  and  the  same  focal 
mechanism,  which  of  course  never  happens.  While  we  have  not  been  able 
to  reproduce  the  high  correlation  levels  found  by  Chael  and  Kromer  (1988) 
(who  analyzed  also  sequence  A  but  not  necessarily  the  same  aftershocks),  we 
have  tried  to  relocate  the  aftershocks  using  a  master-event  technique.  Within 
the  limits  of  the  quality  of  the  relative  arrival  times,  however,  no  significant 
difference  was  found.  This  is  not  surprising  considering  the  fact  that  the  source 
radii  for  the  main  A  and  B  events  have  been  estimated  to  1-2  km  (Hansen 
et  al.,  1989).  Compared  to  the  hypocentral  variations  accepted  by  Boore  and 
Atkinson  (1989),  for  example,  the  data  used  in  this  study  are  therefore  closer 
to  the  ideal  requirement.  An  interesting  and  yet  unresolved  question  in  this 
respect  is,  by  the  way,  what  limits  in  hypocentral  variations  can  be  accepted 
in  studies  of  this  kind? 

The  hedgehog  inversion  approach  that  we  have  used  in  the  present  study 
clearly  requires  high-quality  data  as  well  as  a  careful  data  analysis  procedure. 
The  spectral  smoothing  is  an  essential  part  of  this  procedure.  An  interesting 
option  in  this  inversion  is  that  the  seismic  moments  can  be  (re-)determined 
along  with  the  scaling  parameters,  without  essentially  affecting  the  resolving 
power  with  respect  to  those  parameters.  The  good  correlation  that  we  have 
observed  here  for  both  of  the  earthquake  sequences  reflect  a  high  quality  in 
the  initial  moment  determinations. 

Given  the  determination  of  a  preferred  7  value  (such  as  2).  this  inversion 
determines  simultaneously  a  set  of  6  and  /o  values.  The  shape  and  the  gradi- 
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ents  of  the  contour  levels  in  such  plots  provide  important  information  about 
the  resolving  power,  including  trade-off  effects  between  the  parameters. 

We  have  chosen  in  this  study  to  invert  for,  in  addition  to  6,  the  reference 
event  corner  frequency  /o  rather  than  the  stress  drop.  The  reason  for  this  is 
simply  that  it  is  the  corner  frequency  which  enters  the  spectral  ratio  relation 
(equation  (1)),  and  that  a  jnnection  to  stress  drop  requires  some  additional 
constraint  (model)  for  the  relation  between  corner  frequency  and  stress  drop. 
In  using  corner  frequency,  we  obtain  therefore  primarily  the  stress  drop  ratios, 
as  shown  in  Table  2.1.2.  It  is  important  to  note  here  that  those  values  (as  well 
as  the  corner  frequency  ratios)  are  quite  reliably  determined,  in  contrast  to 
the  absolute  values  of  the  main  event  corner  frequencies.  As  an  example,  the 
latter  values  have  been  used,  in  the  last  column  of  Table  2.1.2,  in  deriving  a 
set  of  stress  drop  values  based  on  the  Brune  source  model. 

These  stress  drop  values  reflect  directly  that  the  corner  frequencies  de¬ 
termined  here  are,  as  mentioned  above,  somewhat  higher  than  those  esti¬ 
mated  earlier,  using  different  source  spectral  approaches  (Hansen  et  al.,  1989; 
Kvamme  and  Hansen,  1991).  While  we  cannot  explain  this  difference  at  this 
stage  other  than  in  general  terms  (lack  of  near-field  observations,  differences 
in  models  and  approaches,  resolution  problems),  there  is  an  important  im¬ 
plication  with  respect  to  the  use  of  the  spectral  ratio  method.  This  is  that 
the  use  of  a  fixed  pre-determined  corner  frequency,  such  as  one  derived  from 
magnitude  and  some  standard  model,  can  easily  bring  us  off  scale  with  respect 
to  the  fit  between  computed  and  observed  spectral  ratios.  In  turn,  this  can 
result  in  biased  6  values,  equivalent  to  traversing  along  the  ridge  crests  in  Fig¬ 
ure  2.1.3.  If  we  replace  the  corner  frequency  with  stress  drop,  the  same  result 
holds  true. 

We  have  found  in  this  study,  when  assuming  an  u2  model,  that  the  stress 
drop  scaling  parameter  6  has  values  of  3.24  and  3.80  for  two  earthquake  se¬ 
quences  located  essentially  within  similar  tectonic  environments.  An  u>3  model 
gives  a  corresponding  difference  with  respect  to  6.  The  potential  implication 
of  this  is  that  we  may  have  to  accept  quite  some  variation  in  the  stress  drop 
parameter,  not  only  between  widely  separated  main  shocks  (Hanks,  1977),  but 
even  between  closely  located  events  and  event  sequences.  Such  observations 
may  serve  to  illustrate  why  the  enduring  discussions  on  the  CSD/ISD  alterna¬ 
tives  (Kanamori  and  Anderson,  1975;  Nuttli,  1983a;  1983b)  and  on  the  possi¬ 
ble  connection  between  stress  drop  and  tectonic  environment  (Sommerville  et 
al.t  1987;  Kanamori,  1988)  still  are  essentially  unresolved.  The  large  scatter  in 
the  data  that  has  nourished  this  discussion  indicates  that  the  intra-earthquake 
variability  within  selected  regions  (such  as  California  and  Fennoscandia)  may 
be  as  large  as  the  variability  between  the  regions.  That  means,  in  turn,  that 
it  may  be  possible  to  use  results  about  earthquake  source  properties  across 
regional  and  tectonic  boundaries  to  a  larger  extent  than  previously  assumed 
(cf.  Campbell,  1989),  and  that  a  large  and  yet  unpredictable  variability  is  still 
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present  in  such  data,  implying  large  uncertainties  with  respect  to  prediction 
of  ground  motion.  As  is  well  known,  such  predictions  are  very  sensitive  to  the 
v  stress  drc  p  parameter. 

A  final  point  here  is  that  the  scaling  relations  derived  in  this  study  should 
not  necessarily  be  assumed  to  be  valid  for  any  magnitude  range  outside  that 
covered  by  the  present  data  (Mi  1.5  to  5).  Our  results  favor  a  single  model 
covering  this  range,  including  also,  as  noted  initially,  the  magnitude  (around  3) 
where  departures  from  a  CSD  situation  often  are  found  to  occur,  for  series  of 
main  earthquakes.  The  results  of  McGarr  (1986)  indicates  that  the  CSD/ISD 
transition  for  earthquake  sequences  occurs,  if  at  all,  at  higher  magnitudes. 
It  seems  therefore  that  scaling  relations  in  general  should  not  without  good 
reasons  be  assumed  to  be  valid  for  data  other  than  those  from  which  the 
relations  have  been  derived. 

Conclusions 

Short-period  and  high-frequency  NORESS  observations  from  two  after¬ 
shock  sequences  offshore  western  Norway,  with  epicentral  distances  of  424  and 
405  km,  respectively,  have  been  used  in  testing  a  new  hedgehog  inversion  pro¬ 
cedure  for  determination  of  source  spectral  scaling  characteristics,  based  on  a 
previously  developed  spectral  ratio  method.  This  inversion  procedure  provides 
insight  into  resolution  properties,  uncertainties,  and  trade-off  considerations 
between  source  parameters. 

The  two  main  events  analyzed  have  Mi  magnitudes  of  about  5  while  the 
smallest  aftershocks  are  about  1.5,  with  seismic  moments  in  the  range  1017  to 
1012  Nm.  The  observed  spectral  ratios  support  w2  source  models  for  both  of 
the  earthquake  sequences,  but  only  if  certain  constraints  are  imposed  on  the 
scaling  parameter  (S)  that  determines  the  stress  drop  sensitivity  with  respect 
to  seismic  moment.  In  that  case,  stress  drop  is  found  to  be  almost  independent 
of  seismic  moment  (6  of  3.24)  for  one  of  the  sequences  and  increasing  with 
moment  (6  of  3.80)  for  the  other.  Without  such  constraints,  which  may  not 
be  justifiable  in  the  case  of  earthquake  sequences,  equally  good  solutions  are 
found  for  w3  models,  in  which  case  6  attains  values  around  5-6,  indicating 
strongly  increasing  stress  drop.  External  constraints  on  the  source  spectral 
slope  and/or  corner  frequency  will  similarly  help  in  resolving  6. 

The  method,  which  also  simultaneously  determines  corner  frequencies  (or 
%  stress  drop  if  adopting  some  model  relating  this  to  corner  frequency)  for  all 

of  the  events,  albeit  with  a  poorer  resolution,  can  either  use  independently 
determined  seismic  moments  or  it  can  determine  these  as  part  of  the  inversion. 
m  Quite  similar  results  are  obtained  in  those  two  cases,  reflecting  essentially  the 

quality  of  the  initial  moment  determinations.  The  observed  variations  in  stress 
drop  sensitivity  between  two  earthquake  sequences,  only  84  km  apart  and 
occurring  in  similar  tectonic  environments,  are  consistent  with  the  considerable 
variations  that  have  been  observed  on  a  larger  scale.  Caution  is  therefore 
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needed  if  the  scaling  relations  derived  in  this  study  are  applied  to  other  events 
(larger  ones  and/or  series  of  main  events),  or  to  different  seismic  regions. 


H.  B  ungum 
A.  Alsaker 
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Event 

Year 

Date 

Hour 

Data 

Ml 

No 

Mo 

(Nm) 

No 

Ao 

1986 

Feb 

05 

17.53.35 

hf-sp 

4.8 

1 

6.5  • 

1015 

2 

Ai 

20.23.15 

hf-sp 

2.2 

1 

4.7 

1012 

1 

a2 

23.35.42 

hf-sp 

1.8 

1 

2.5  • 

1012 

1 

a3 

Feb 

06 

06.19.51 

hf-sp 

1.8 

1 

1.4 

1012 

1 

Aa 

Feb 

10 

12.31.36 

hf-sp 

1.4 

1 

1.5  ■ 

1012 

1 

Ai 

Feb 

13 

13.39.01 

sp 

2.0 

1 

3.1 

1012 

1 

A6 

Feb 

13 

19.03.50 

sp 

2.0 

1 

3.5- 

1012 

1 

Bo 

1989 

Jan 

23 

14.06.29 

hf-sp 

4.9 

6 

4.3- 

1016 

6 

B\ 

14.35.21 

hf-sp 

1.6 

1 

2.6- 

1012 

1 

b2 

15.57.27 

hf-sp 

1.9 

5 

2.8- 

1012 

4 

Bz 

16.40.31 

hf-sp 

2.8 

4 

2.6- 

1013 

4 

Ba 

21.38.05 

hf-sp 

2.5 

6 

1.3  • 

1013 

6 

Bs 

Jan 

25 

23.07.36 

hf-sp 

1.5 

4 

1.2- 

1012 

3 

Be 

Jan 

26 

10.01.43 

hf-sp 

1.6 

4 

1.8- 

1012 

2 

Table  2.2.1.  Events  analyzed  in  this  study,  with  event  code,  year,  date, 
origin  time,  NORESS  data  channels,  local  magnitude  (Ml),  number  of  records 
used  in  estimating  Ml,  predetermined  seismic  moments  (Nm),  and  number  of 
records  used  in  estimating  Mq. 


Event 

Mo 

(Nm) 

fo  (Hz) 

oy/<r,- 

1 7  (MPa) 

Ao 

(1.8 

1016) 

2.19 

(1.00) 

10.4 

Ai 

2.1 

1013 

18.34 

0.60 

6.3 

a2 

8.1 

1012 

24.64 

0.56 

5.9 

A3 

3.4 

1012 

32.50 

0.53 

5.5 

Aa 

4.2 

1012 

30.30 

0.54 

5.6 

■d.5 

1.2 

1013 

22.05 

0.58 

6.0 

>46 

1.2 

1013 

21.77 

0.58 

6.1 

Bo 

(7.0 

1016) 

1.53 

(1.00) 

23.4 

Bi 

4.9 

1012 

18.96 

0.13 

3.1 

b2 

3.9 

1012 

20.03 

0.13 

3.0 

Bz 

8.2 

1013 

9.00 

0.24 

5.7 

Ba 

3.9 

1013 

10.99 

0.21 

4.8 

Bz 

1.3 

1012 

26.78 

0.10 

2.4 

Bo 

1.8 

1012 

24.76 

0.11 

2.5 

Table  2.1.2.  Results  from  the  present  source  scaling  inversion  analysis,  with 
event  code,  seismic  moment  (not  determined  for  Ao  and  B0),  corner  frequency, 
stress  drop  ratios,  and  stress  drops  when  using  equation  (5)  for  determining 
the  Brune  (1970,  1971)  stress  drop  for  the  reference  event  (1  MPa  is  10  bars). 
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Fig.  2.1.1.  Observed  spectral  ratios  for  event  sequence  A  (fully  drawn  lines) 
together  with  theoretical  ones  (dashed  lines)  for  the  set  of  parameters  (see 
top  label)  that  gives  the  best  fit,  including  also  an  adjustment  of  the  moment 
ratios.  Numerical  values  for  this  solution  are  given  in  Table  2.1.2. 


Frequency  (Hz) 

Fig.  2.1.2.  Spectral  ratios  as  in  Figure  2.1.1,  but  for  earthquake  sequence  B. 
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Scaling  parameter  <5  Scaling  parameter 


B  Inversion  including  moment  ratios 


*o 


4.4 

4.0 

3.6 

3.2 

2.8 


4.0 

3.9 

3.8 

3.7 

3.6 


1.4  1.6 

Corner  frequency  (Hz) 


1.40  1.50  1.60 

Corner  frequency  (Hz) 


Fig.  2.1.3.  Contour  diagrams  of  root-mean-square  (RMS)  deviations  between 
computed  and  observed  spectral  ratios  for  event  sequence  B,  as  a  function 
of  reference  event  corner  frequency  /o,r  and  the  scaling  parameter  S.  An  u>2 
source  model  (7  of  2)  is  assumed  here,  and  the  four  diagrams  cover  the  first  four 
iterations,  where  moment  ratios  are  adjusted  simultaneously.  The  numbers  1 
to  4  in  each  frame  indicate  iteration  number  and  also  the  locations  of  the 
minima. 
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B  Old  moments  7=3  B  new  moments  7=3 


Fig.  2.1.4.  Contour  diagrams  as  for  Figure  3.1.2,  but  for  an  u3  source  model 
(7  of  3).  The  parameter  ranges  covered  here  are  twice  as  large  as  the  earlier 
ranges,  which  are  indicated  by  the  boxes  in  the  lower-left  corners  of  each  frame. 


Corner  frequency  (Hz)  Corner  frequency  (Hz) 


Fig.  2.1.5.  Trajectories  of  source  scaling  inversion  solutions  when  the  source 
spectral  shape  parameter  7  is  allowed  to  vary  from  1.5  to  3.0,  in  steps  of  0.1. 
The  values  in  the  parentheses  are  7  and  resulting  RMS  value  corresponding 
to  that  solution,  and  the  two  dashed  lines  indicate  the  initial  range  considered 
for  the  parameter  6  (3  is  CSD,  4  is  ISD).  The  left  frame  is  for  earthquake 
sequence  A  and  the  right  one  is  for  B. 
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2.2  P-wave  focusing  effects  at  NORSAR  for  Novaya  Zemlya 
explosions 

Introduction 

It  is  now  well  established  that  significant  improvements  in  P-wave  de¬ 
tectability  of  weak  seismic  events  can  be  achieved  by  taking  advantage  of  re¬ 
ceiver  focusing  effects  of  P-waves.  Such  effects  appear  to  originate  from  upper 
mantle  heterogeneities  underneath  the  receiver,  and  are  strongly  dependent 
upon  incoming  azimuth  and  angle  of  incidence  of  the  wavefront.  Thus,  at 
any  given  site,  P-wave  focusing  is  expected  to  be  strong  for  certain  epicentral 
regions,  and  weak  for  others,  causing  highly  variable  sensitivity  for  P-wave 
detection  of  signals  from  different  source  regions. 

One  of  the  most  well-known  examples  of  P-wave  focusing  is  the  remark¬ 
able  detection  performance  of  the  NO  RES  S  array  for  Semipalatinsk  explosions 
(Ringdal,  1990).  NORESS  is  located  at  a  site  (06C02)  within  the  large  NOR¬ 
SAR  array  that  has  especially  favorable  focusing  effects  for  P-waves  from  the 
Semipalatinsk  region,  and  detectability  has  been  further  amplified  through 
the  array  gain  of  14-16  dB  in  the  frequency  band  (2-4  Hz)  of  highest  signal- 
to-noise  ratio.  On  the  other  hand,  for  many  other  source  areas,  the  focusing 
effects  at  the  NORESS  site  are  far  less  favorable,  and  the  array  detection 
performance  is  correspondingly  lower. 

In  this  paper  we  study  the  P-wave  focusing  effects  across  NORSAR  for 
Novaya  Zemlya  explosions,  and  give  some  preliminary  comments  on  the  results 
achieved  so  far. 

NORSAR  P-wave  recordings  of  Novaya  Zemlya  explosions 

A  number  of  underground  nuclear  explosions  have  been  conducted  at  No¬ 
vaya  Zemlya  during  the  20  years  of  NORSAR  operation.  Because  of  the  prox¬ 
imity  of  NORSAR  to  this  test  site  (20  degrees)  and  the  excellent  P-wave 
propagation  characteristics,  the  recorded  signals  have  usually  been  so  strong 
that  they  have  exceeded  the  dynamic  range  of  the  NORSAR  digital  recording 
system.  In  fact,  even  the  smallest  explosion  reported  by  Lilwall  and  Marshall 
(1986),  (9  October  1977,  mj,  =  4.51)  was  in  this  category.  Fig.  2.2.1  shows 
NORSAR  recordings  (7  SP  seismometers,  1  per  subarray)  for  this  event.  It 
can  be  noted  that  two  of  the  data  channels  are  clipped  (02C01  and  03C01). 
In  the  latter  case  (03C01)  the  amount  of  clipping  is  very  significant. 

This  example  shows  that  for  NORSAR  seismometers  located  at  sites  with 
favorable  receiver  focusing  effects,  even  very  small  explosions  (mj,  >  4.0)  at 
Novaya  Zemlya  produce  high  amplitude  recordings.  While  this  is  a  disad¬ 
vantage  for  studying  focusing  effects  and  other  waveform  features,  it  clearly 
implies  that  there  is  a  very  significant  potential  for  exploiting  these  effects  in 
achieving  detectability  at  low  magnitude  levels. 
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Among  the  presumed  explosions  at  Novaya  Zemlya  recorded  by  NORSAR, 
we  have  found  only  two  events  that  are  suitable  for  an  amplitude  pattern  study 
(Events  2  and  3  of  Table  2.2.1). 

The  location  estimates  in  Table  2.2.1  of  these  two  events  (2  and  3)  are 
approximate,  but  both  explosions  are  from  the  Northern  test  site.  The  mj, 
values  have  been  estimated  as  world-wide  equivalent,  by  comparing  NORSAR 
recordings  to  Event  1  (9  October  1977). 

Figures  2.2.2  and  2.2.3  show  NORSAR  P-wave  recordings  of  events  2  and 
3,  using  the  same  instrument  selection  as  for  Event  1  (Figure  2.2.1).  We  note 
that  on  Figure  2.2.2,  the  background  (microseismic)  noise  level  is  very  high, 
but  the  signal  is  still  clearly  visible  on  some  traces  (notably  03C01).  On  Figure 
2.2.3,  the  general  amplitude  pattern  is  similar  to  the  other  two  events,  again 
with  02C01  and  03C01  showing  the  strongest  signals. 

NORSAR  amplitude  patterns 

In  practice,  detectability  is  determined  by  the  signal-to-noise  ratio  in  the 
“best”  frequency  band.  For  NORSAR  recordings  of  Novaya  Zemlya  explo¬ 
sions,  the  2.5-4.5  Hz  band  is  near  optimum,  and  Figure  2.2.4  shows  the  data 
from  Figure  2.2.2  after  applying  this  filter.  The  amplitude  pattern  in  Figure 
2.2.4a  (given  by  the  scaling  factors  to  the  left  of  each  trace)  shows  an  order 
of  magnitude  variation  across  the  60  km  aperture  array.  From  Figure  2.2.4b 
it  can  be  seen  that  the  noise  levels  at  the  different  sites  are  similar,  thus  the 
signal  amplitude  pattern  reflects  true  differences  in  P-wave  detectability. 

There  are  also  significant  amplitude  variations  within  NORSAR  subarrays 
(typical  diameter  8-10  km),  although  naturally  not  as  great  as  across  the 
entire  NORSAR  array.  Figures  2.2.5  and  2.2.6  show  as  examples  unfiltered 
data  for  event  3  within  the  two  best  subarrays  (02C  and  03C).  Note  that  the 
pattern  is  generally  stable  within  subarray  02C,  while  at  03C  the  previously 
selected  site  (03C01)  stands  out  as  having  far  stronger  signal  than  the  others. 

Figure  2.2.7  gives  an  overall  view  of  NORSAR  P-wave  amplitude  levels  for 
Novaya  Zemlya  explosions  as  a  function  of  seismometer  location.  It  is  clear 
that  there  axe  many  sites  with  very  favorable  focusing  effects,  in  particular 
O3C0I.  Compared  to  the  NORESS  site  (06C02),  there  is  an  order  of  magnitude 
difference. 

Detectability  considerations 

Across  the  large  NORSAR  array,  the  coherency  of  the  signals  from  Novaya 
Zemlya  is  far  too  low  to  allow  any  meaningful  improvement  in  signal-to-noise 
ratio  through  standard  array  processing.  In  fact,  the  NORSAR  array  beam 
has  a  lower  SNR  than  the  best  single  sensors  for  this  region. 

Thus,  when  considering  seismic  events  in  the  regional  distance  range,  the 


22 


main  benefit  of  having  full  NORSAR  data  available  is  to  obtain  location  esti¬ 
mates  for  small  events,  and,  very  importantly,  to  identify  sites  with  particu¬ 
larly  strong  signals  for  a  given  region  of  interest.  It  is  easy  (taking  appropriate 
reservations  in  view  of  the  limited  data  available)  to  make  first-order  estimates 
of  what  could  be  achieved  in  this  way.  E.g.,  by  simple  SNR  scaling,  we  have 
found  that  the  seismometer  03C01  would  record  detectable  signals  for  Novaya 
Zemlya  explosions  at  the  northern  test  site  of  m&  =  2.0-2. 5  during  normal 
noise  conditions.  With  a  small  array  of  6-10  instruments  at  that  site  (diam¬ 
eter  ca  1000  m),  it  is  projected  that  a  further  improvement  of  0.5  m*  units 
could  be  achieved.  This  is  based  on  NORESS  noise  structure  experience,  and 
would  imply  a  threshold  below  m*,  =  2.0  for  such  a  small  array.  Again,  it  is 
emphasized  that  further  studies  are  necessary  before  the  potential  for  such  low 
level  detection  can  be  confidently  established. 


F.  Ringdal 
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Event  1  9  Oct  1977  73.414°N  54.935°E  m6  =  4.51  O.T.  10.59.58.1 

Event  2  15  Nov  1978  73  °N  55  °E  mb  =  3.6  O.T.  08.30.00.0 

Event  3  26  Aug  1984  73  °N  55  °E  mb  =  3.8  O.T.  03.30.00.0 


Table  2.2.1.  Presumed  explosions  at  Novaya  Zemlya  analyzed  in  this  study. 
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Fig.  2.2.1.  NORSAR  unfiltered  recordings  (one  seismometer  for  each  of  7 
subarrays)  for  the  Novaya  Zemlya  underground  nuclear  explosion  of  9  October 
1977  (Event  1  in  Table  2.2.1).  The  scaling  factors  are  shown  to  the  left  of  each 
trace.  Note  that  even  though  this  explosion  was  very  small  (m*  «  4.5),  the 
recorded  signal  was  strong  enough  to  saturate  the  NORSAR  digital  recording 
system  at  sites  02C01  and  03C01. 
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Fig.  2.2.2,  Same  as  Fig.  2.2.1,  but  for  Event  2  (m^  «  3.6)  in  l^ble  2.2.1. 
Note  that  in  spite  of  very  strong  background  noise,  the  signal  is  clearly  visible 
at  sites  02C01  and  03C01,  even  on  the  unfiltered  recordings. 
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Pig.  2.2.3.  Same  as  Fig.  2.2.1,  but  for  Event  3  (mj  ss  3.8 )  in  Table  2.2.1. 
Note  that  the  background  noise  in  this  case  is  much  lower  than  on  Fig.  2.2.2. 
Again  sites  02C  and  03C  have  the  strongest  signals. 
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Fig.  2.2.4a.  Filtered  recordings  (2.5-4. 5  Hz)  of  the  data  shown  in  Fig.  2.2.2. 
Note  that  the  signal  can  be  seen  on  all  traces,  but  with  maximum  amplitudes 
varying  by  more  than  an  order  of  magnitude  (see  scaling  factors  to  the  left  of 
each  trace). 
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Fig.  2.2.4b.  Noise  data  (2.5-4.5  Hz)  preceding  the  signal  shown  in  Fig.  2.2.4a. 
All  traces  have  the  same  scale.  Note  that  the  noise  levels  at  different  sites  are 
similar,  thus  the  amplitude  variation  in  Fig.  2.2.4a  reflects  true  differences  in 
P-wave  detectability  of  Novaya  Zemlya  explosions. 
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Fig.  2.2.5.  P-wave  amplitude  variation  within  subarray  03C  for  Event  3 
(unfiltered  data).  There  is  a  factor  of  6  in  variation  from  the  weakest  to  the 
strongest  signal. 
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Fig,  2.?.0.  P-wave  amplitude  variation  within  subarray  02C  for  Event  3 
(unfiltered  data).  Note  that  both  the  amplitudes  and  waveform  shapes  are 
similar  within  this  subarray,  in  contrast  to  subarray  03C  (Fig.  2.2.5). 


NORSAR  amplitude  pattern  -  Novaya  Zemlya 
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Fig.  2.2.7.  P-wave  magnitude  bias  relative  to  average  NORSAR  m 4  for  indi¬ 
vidual  sites  within  the  array.  Symbol  sizes  are  proportional  to  nu  bias,  with 
plusses  indicating  positive  bias  values.  Note  the  high  positive  bias  for  all  sites 
within  subarray  02C,  and  the  especially  high  bias  value  at  03C01.  The  range 
of  the  bias  values  is  from  +0.9  (03C01)  to  -0.3  (01B05). 
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2.3  Wavefield  extrapolation  of  array  data 

A  major  consideration  in  selecting  the  NORESS  array  site  has  been  the  high 
amplification  of  short-period  P  signals  from  directions  to  the  East.  From  the 
analysis  of  data  from  large  arrays  it  is  well-known  that  the  amplitude  anoma¬ 
lies  are  associated  with  significant  time  residuals  (Berteussen,  1976).  It  is 
desirable  and  possible  to  calibrate  the  times  and  amplitudes  of  specific  phases 
at  NORESS  and  similar  arrays,  but  it  would  be  an  advantage  both  from  a 
scientific  and  from  a  practical  point  of  view  to  base  the  calibration  and  inter¬ 
polation  of  seismic  anomalies  on  a  physically  acceptable  model  of  subsurface 
structure.  As  a  first  step  toward  that  goal,  we  describe  an  experiment  to 
extrapolate  the  observed  wavefield  at  the  surface  down  onto  the  subsurface. 
Downward  extrapolation  is  common  in  processing  seismic  reflection  data  where 
the  procedure  is  part  of  a  migration  process.  The  experiment  we  describe  here 
differs  from  the  usual  reflection  methods  in  two  important  aspects:  (1)  The 
time  and  amplitude  anomalies  are  the  result  of  forward  scattering,  rather  than 
backward  scattering  as  in  reflection  seismology.  Downward  extrapolation  of 
the  forward  scattered  wavefield  has  sometimes  been  called  holography  (Troit- 
skiy  et  al,  1981);  (2)  The  wavefield  at  the  Earth’s  surface  is  sampled  on  a  more 
or  less  irregular  grid  (c.f.  the  NORSAR  array  configuration).  This  precludes 
the  effective  use  of  Fourier  transform  methods.  However,  the  elastodynamic 
representation  theorem  can  be  used  to  formulate  wavefield  extrapolation  as  an 
inverse  problem  independent  of  how  the  field  is  sampled.  A  range  of  methods 
is  then  available  to  solve  this  problem. 

The  inverse  problem 

Let  the  objective  of  downward  extrapolation  be  to  determine  the  field  on 
subsurface  5.  The  field  at  the  Earth’s  surface  in  x  is  given  by  (c.f.  Doornbos, 
1988): 

« «(£»0  =  Js{T'jk(i,x,t)  *  u j(£,t)  -  G’(£,z,t)  *  Tjk(t,t)}nk  dS  (1) 

Here  £  is  a  point  on  5  with  upward  normal  n,  the  displacement  and  stress 
fields  are  denoted  by  uj  and  rjk,  respectively,  and  the  Green’s  function  and 
its  associated  stress  tensor  by  G|  and  T'-k. 

Although  it  is  customary  to  simplify  equation  (1)  to  a  Kirchhoff  integral 
under  the  assumption  of  acoustic  waves,  we  use  here  a  less  restrictive  far  field 
approximation: 


«;(£,<)  =  A»(£)s;(£)/(t-T,„) 

G)(i,x,t)  =  Ag(£;x)s,(x;£)  s>(£;z)  6(t  -Tq)  (2) 

where  T,n  =  T(£),  Ta  =  T(£;  x),  and  the  unit  vectors  s  indicate  the  direction 
of  particle  motion  in  £  and  x.  Thus  if  v  is  a  unit  wave  motion  vector,  then 
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s-u—  1  for  P,  and  s-u  =  0  for  S  waves.  Substituting  the  asymptotic  forms 
(2),  assuming  the  amplitude  factors  Ain  and  Ac  to  be  real,  and  disregarding 
time  and  amplitude  anomalies,  reduces  equation  (1)  to  a  representation  of  the 
reference  field: 

a°(x,  t)=  f  S(x- 0 B°(£;  x)  f(t  -  T°)  dS  (3) 

Js  ~ 

where 

T°  =  Tin(t)  +  TG(i;x)  (4) 

X)  =  Atn(0  Ag(£;,  X)  R(, f;  x)  (5) 

The  factor  R(£;x)  depends  on  the  incident  and  scattered  wave  modes.  We 
will  be  concerned  with  P  to  P  scattering.  Then  if  subsurface  5  is  embedded 
in  an  isotropic  region  with  Lame  constants  A  and  n  and  P  velocity  a: 

R  =  —  {A  -  2n(v  •  u')}  {(«/  •  n)  -  (u  ■  n)}  (6) 

a 

Based  on  equation  (3),  the  actually  observed  wavefield  u(x,t)  will  now  be 
written  as  a  perturbation  of  an  initial  estimate  at  subsurface  S: 

u(x,  t)  =  [  s  (1  +  6B)  B  f(t  -  T  -  6T)  dS  (7) 

Js 

Here  the  amplitude  and  time  perturbations  SB  and  ST  are  unknown,  but 
B,  T  and  f(t )  are  assumed  known.  Note  that  B  and  T  need  not  be  equal  to 
B°  and  T°  of  equation  (3).  In  general,  B  and  T  should  be  chosen  such  that  the 
perturbations  SB  and  ST  are  relatively  small.  In  practice,  values  for  B  and 
T  can  be  obtained  from  the  results  of  downward  extrapolating  the  wavefield 
to  a  level  slightly  above  S.  The  waveform  f(t)  is  obtained  from  the  array 
beam,  assuming  that  the  beam  provides  an  estimate  of  the  reference  signal 
u°(t)  and  noting  that  u°(t)  ~  f(t  -  T°).  The  representation  (7)  is  physically 
meaningful  if  the  perturbations  are  associated  with  topography  of  5  and/or 
laterally  varying  structure  below  5  (provided  this  structure  does  not  produce 
caustics  below  5). 

Using  equation  (7)  and  the  approximation  for  the  modified  reference  field: 

=  J  sB  j(t  -  t)  dS  (8) 

We  can  pose  a  linearized  inverse  problem  in  3  different  forms: 

(a)  Bom  approximation  in  time  domain: 

6u(x,t)  =  u(x,  t)  -  u(x,  t) 

~  /  s  B{SB  f(t  —  T)  -  ST  f(t  -  T)  dS  (9) 

Js 
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(b)  Bom  approximation  in  frequency  domain: 


ou(x,  w)  ~  —iw  F(w)  J  s  B  e,wT(6B  +  iwST)  dS 


(10) 


(c)  Rytov  approximation: 

. U(x,w )  F(w) 


lnw<rh  -  ~iwW(\  f 1 B  *iwT(<6B  +  iw6T)  ds  (“) 

U(x,w)  U(x,w)Js 


Synthetic  data 

We  have  synthesized  the  surface  displacement  u(x,  t)  from  prescribed  time 
and  amplitude  anomalies  at  depth.  The  spatial  sampling  of  u(x,  t )  was  chosen 
so  as  to  simulate  the  recordings  of  132  short-period  channels  in  the  NORSAR 
array.  The  array  beam  signal  from  a  nuclear  explosion  in  E.  Kazakh  was 
used  to  get  the  basic  waveform  /(t)  and  the  direction  of  wave  propagation 
v.  Downward  extrapolation  was  then  done  in  3  different  ways,  by  inverting 
equations  (9),  (10)  and  (11),  respectively.  The  time  domain  method  (a)  used 
the  first  3  seconds  of  the  low-pass  filtered  P  pulse  resampled  at  5  Hz.  The 
frequency  domain  methods  (b)  and  (c)  used  6  frequencies  of  the  same  pulse, 
in  the  band  0.6  -  2.1  Hz. 

Figure  2.3.1  pictures  an  anomaly  at  100  km  depth  used  in  one  of  the 
experiments.  In  this  example  the  anomaly  equals  a  time  delay  with  maximum 
0.7  seconds.  The  resulting  time  and  amplitude  patterns  at  the  NORSAR  array 
site  are  shown  in  Figures  2.3.2  and  2.3.3.  The  time  residuals  are  reduced  in 
relation  to  the  pattern  at  100  km  depth,  but  the  amplitudes  fluctuate  by 
about  a  factor  of  2.  Inversion  of  these  data  was  started  at  20  km  depth,  and 
continued  downward  with  a  5  km  depth  step.  For  the  purpose  of  inversion, 
the  subsurface  area  at  each  depth  was  divided  in  cells,  the  size  of  which  is 
increasing  with  distance  from  the  central  ray.  The  reconstructed  pattern  at 
100  km  depth  is  shown  in  Figure  2.3.4  and  2.3.5,  where  we  have  used  the 
Rytov  method  (c).  There  are  outliers  along  the  boundaries  of  the  grid  where 
data  coverage  is  poor.  There  is  also  a  time  offset  of  about  0.2  s.  However, 
the  principal  features  of  the  time  residual  pattern  are  well  recovered.  On  the 
other  hand,  the  amplitudes  fluctuate  by  about  a  factor  of  2  in  contrast  to  the 
starting  model;  we  interpret  this  to  mean  that  the  surface  data  are  primarily 
constraining  the  time  perturbations  at  depths. 

Explosion  P  wave  from  E.  Kazakh 

Figure  2.3.6  shows  the  P  wave  on  the  C-ring  of  NORSAR,  from  a  nuclear 
explosion  in  E.  Kazakh.  The  time  residuals  vary  by  more  than  0.5  s,  and 
the  amplitudes  vary  by  about  a  factor  of  5  (Figures  2.3.7  and  2.3.8).  The 
data  from  the  complete  array  were  used  to  reconstruct  the  residual  time  and 
amplitude  pattern  at  depth.  Results  for  100  km  depth  are  shown  in  Figures 
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2.3.9  and  2.3.10.  From  the  results  for  synthetic  data  we  expect  solutions 
for  the  residual  times  to  be  more  reliable  than  those  for  the  amplitudes.  The 
variation  of  the  residual  times  is  about  twice  as  large  as  at  the  Earth’s  surface. 
Disregarding  the  poorly  constrained  outliers  along  the  boundaries  of  the  grid, 
the  results  suggest  an  extensive  deep-seated  region  of  low  velocity  to  the  NE 
from  NORSAR  (note  that  the  center  of  the  maps  2.3.9  and  2.3.10  is  already 
displaced  to  the  NE  from  the  array  center).  In  this  connection  it  should  also 
be  noted  that  the  optimum  depth  of  subsurface  S  is  actually  another  unknown 
of  the  inverse  problem.  Several  error  criteria  may  be  used  to  determine  this 
depth.  In  Figure  2.3.11  we  have  plotted  the  RMS  error  of  time  residuals 
left  unexplained  by  the  subsurface  solution,  in  per  cent  of  the  observed  RMS 
residual.  Inversion  of  the  synthetic  data  reduces  the  RMS  error  to  very  small 
values  for  depths  larger  than  about  40  km.  Inversion  of  the  explosion  data 
leaves  a  more  significant  part  of  the  residuals  unexplained,  but  there  is  a  broad 
minimum  in  the  range  100-130  km.  It  is  of  course  desirable  to  delineate  the 
structure  in  more  detail;  inversion  of  the  data  from  many  more  events  is  needed 
to  provide  this  detail. 

E.  0degaard,  Univ.  of  Oslo 
D.J.  Doornbos,  Univ.  of  Oslo 
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Fig.  2.3.1.  Contour  map  of  travel  time  delays  at  100  km  depth.  The  center 
of  the  map  is  at  60.992°N,  12.134°E  (NE  from  the  NORSAR  array  center). 
The  incident  ray  slowness  and  propagation  azimuth  are  0.076  s/km  and  255 
deg. 
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Fig.  2.3.2.  Travel  time  residuals  at  the  Earth’s  surface,  caused  by  the 
anomaly  of  Fig.  2.3.1.  The  NORSAR  subarray  centers  are  indicated  by  stars. 
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100  km 


Fig.  2.3.3.  Amplitude  pattern  at  the  Earth’s  surface,  caused  by  the  anomaly 
of  Fig.  2.3.1.  Amplitude  factors  relative  to  the  reference  value  1. 
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100  km 


Fig.  2.3.4.  Reconstructed  time  residuals  at  100  km  depth,  the  map  covers 
the  same  area  as  Fig.  2.3.1. 
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Fig.  2.3.6.  NORSAR  subarray  beams  at  C01-C14.  Scaling  factors  to  the 
left.  The  P  wave  is  from  a  nuclear  explosion  in  E.  Kazakh,  in  1976.  Start  time 
of  the  records  underneath. 
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Fig.  2.3.7.  Travel  time  residuals  at  NORSAR,  of  P  wave  from  E.  Kazakh 
nuclear  explosion  (Fig.  2.3.6). 
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Fig.  2.3.8.  Amplitude  factors  at  NORSAR,  of  P  wave  from  E.  Kazakh  nuclear 
explosion  (Fig.  2.3.6). 
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100  km 


Fig.  2.3.9.  Reconstructed  time  residuals  at  100  km  depth.  The  map  center 
is  at  60.992°N,  12  1°4°E  (as  in  Fig.  2.3.1. 
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100  km 


Fig.  2.3.10.  Reconstructed  amplitude  pattern  at  100  km  depth.  The  map 
covers  the  same  area  as  Fig.  2.3.9. 
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Surface  time  residual  error  (21 


Fig.  2.3.11.  RMS  residual  time  error  after  correction  for  reconstructed  time 
and  amplitude  models  at  depth.  Values  in  per  cent  of  RMS  residual  with 
respect  to  the  reference  model.  A,  B  and  C:  Downward  extrapolation  of  syn¬ 
thetic  data.  D:  Downward  extrapolation  of  P  wave  from  nuclear  explosion  in 
E.  Kazakh  (Fig.  2.3.6). 
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2.4  The  seismicity  of  Norway  and  surrounding  areas 
Introduction 

Norway,  the  epicontinental  North  Sea  and  Barents  Sea,  and  the  Norwe¬ 
gian  continental  margin  are  characterized  by  low-to-intermediate  seismicity 
(Kvamme  and  Hansen,  1989,  Hansen  et  al,  1989;  Bungum  et  al ,  1991).  The 
geology  of  the  areas  south  of  62°N,  and  in  particular  the  northern  North  Sea, 
has  been  widely  investigated  as  a  result  of  hydrocarbon  exploration.  The  areas 
north  of  62°N  have  been  extensively  investigated  as  well,  although  not  in  as 
much  detail  as  further  south.  The  regional  geology  and  tectonic  setting  in  all 
of  these  areas  are  therefore  now  relatively  well  understood  (e.g.  Spencer  et  al, 
1984,  1986;  Glennie,  1986). 

Norway  and  its  surrounding  areas,  in  particular  the  northern  continental 
margin,  have  until  recently  been  relatively  poorly  covered  in  terms  of  seis¬ 
mic  instrumentation  (Husebye  et  al,  1975,  1978;  Bungum  and  Fyen,  1980). 
In  the  last  decade  the  situation  has  been  much  improved,  however,  through 
the  installation  of  new  regional  and  local  networks  (Bungum  et  al,  1986). 
This  increased  number  of  stations  has  resulted  in  more  detailed  delineation 
of  seismicity  patterns  as  well  as  improved  focal  mechanism  solutions.  Such 
information  can  in  turn  be  correlated  with  geologic  information. 

One  of  the  incentives  for  the  improvements  in  our  knowledge  of  the  seis- 
motectonics  of  Norway  and  the  Norwegian  continental  shelf  has  been  the  need 
for  better  and  more  reliable  earthquake  hazard  analyses  for  offshore  industrial 
installations.  These  goals  have  been  pursued  through  the  installation  of  local 
microearthquake  networks  (Bungum  et  al,  1986;  Kvamme  and  Hansen,  1989; 
Havskov  et  al,  1991),  and  through  research  studies  specifically  aimed  at  earth¬ 
quake  hazard  and  loading  (Bungum  and  Seines,  1988).  For  offshore  Norway, 
known  in  great  detail  geologically,  such  networks  and  seismicity  studies  are 
moreover  contributing  significantly  to  our  understanding  of  deep  structures 
and  large-scale  features  such  as  the  Viking  Graben. 

Another  important  incentive  for  investigating  the  natural  seismicity  of  Nor¬ 
way  and  surrounding  areas  is  that  the  nuclear  explosion  monitoring  programs 
need  the  best  possible  background  information  about  the  seismicity,  not  only 
in  the  areas  where  such  explosions  could  take  place,  but  also  in  the  areas  within 
local  and  regional  distance  from  the  monitoring  stations  (Richards,  1988).  In 
the  new  Intelligent  Monitoring  System  (Bache  et  al,  1990;  Bratt  et  al,  1990) 
now  installed  at  NORSAR  such  background  seismicity  can  be  displayed  as 
overlays.  For  this  purpose,  three  different  catalogues  have  been  prepared,  cov¬ 
ering  the  three  time  periods  1880-1954,  1955-1979  and  1980-1989,  shown  in 
Figures  2.4.1-2.4.3,  respectively.  The  purpose  of  this  presentation  is  to  docu¬ 
ment  these  catalogues,  and  to  outline  briefly  their  seismotectonic  context. 
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Earthquake  data 

Even  though  the  seismicity  of  Norway  is  moderate,  historical  records  show 
(Bungum  and  Seines,  1988)  that  earthquakes  strong  enough  to  be  felt  clearly 
over  several  hundred  kilometers  have  occurred  repeatedly,  particularly  along 
the  west  coast.  Systematic  collection  of  macroseismic  questionnaires  was  ini¬ 
tiated  in  the  1880’s  at  the  Bergen  Museum,  and  such  data  have  been  collected 
more  or  less  uninterrupted  for  these  areas  ever  since  (Muir  Wood  et  al,  1388). 
This  has  provided,  in  spite  of  a  sparse  population,  a  good  basis  for  a  major 
reassessment  performed  recently  of  the  historical  seismicity  in  Norway  and 
surrounding  areas  (Ambraseys,  1985;  Muir  Wood  and  Woo,  1987).  This  re¬ 
assessment  has  primarily  included  macroseismic  intensity  maps  (“felt  areas”), 
evaluated  in  a  consistent  way  for  this  time  period  by  going  back  to  the  original 
data  for  all  of  the  earthquakes.  From  this  felt  area  information  similarly  con¬ 
sistent  magnitudes  (Ms)  have  been  evaluated,  which  in  turn  have  been  very 
important  for  assessments  of  earthquake  hazard  in  these  areas  (Bungum  and 
Seines,  1988). 

This  historical  earthquake  catalogue  has  been  updated  to  the  end  of  1989 
by  including  all  events  that  could  be  evaluated  on  the  basis  of  felt  effects.  In 
addition,  information  on  the  regional  earthquake  distribution  in  Fennoscandia 
has  now  been  available  for  several  tens  of  years  based  on  data  from  global  net¬ 
works  of  seismological  stations.  The  earlier  global  networks,  however,  provided 
rather  poor  coverage  in  this  area  where  the  seismicity  is  low- to- intermediate 
(Bungum  et  al,  1986;  Bungum  et  al,  1991).  Only  for  the  time  period  since 
1963  do  these  networks  become  important  for  this  region,  and  then  first  of  all 
for  the  area  «  north  of  70°N  (Sykes,  1967).  After  about  1980,  microearthquake 
networks  were  also  installed  in  northwestern  Europe  (Bungum  et  al,  1986; 
Havskov  et  al,  1991),  and  major  improvements  were  again  possible. 

The  historical  data  are  limited  by  the  shortcomings  of  any  macroseismic 
catalogue  in  less  populated  areas,  and  in  particular  for  offshore  areas.  Similar 
limitations  in  the  global  network  solutions  occur  because  very  few  earthquakes 
have  been  large  enough  to  be  detected  by  a  sufficient  number  of  stations.  These 
limitations  have  been  overcome  in  part  in  the  present  work  by  collecting  a 
very  large  number  of  epicentral  solutions  from  a  variety  of  reporting  agencies, 
and  then  carefully  screening  this  data  base  so  as  to  increase  resolution  while 
maintaining  accuracy  (Bungum  and  Seines,  1988;  Bungum  et  al,  1991).  One 
of  the  most  robust  sorting  criterion  in  this  respect  was  found  to  be  the  number 
of  stations  and/or  readings  used  in  computing  the  epicenters. 

Other  sorting  criteria  included  reporting  agency  (not  equally  reliable), 
magnitude,  macroseismic  information,  various  precision  estimates  (including 
consistency  between  different  reports  for  the  same  event),  in  addition  to  num¬ 
ber  of  stations  and  readings  used  in  the  solutions.  A  time  dependent  combi¬ 
nation  of  these  criteria  has  been  used  when  selecting  events  for  the  seismicity 
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maps  in  Figures  2.4.1-2.4.3.  Those  maps  include,  in  addition  to  earlier  data, 
very  significant  contributions  resulting  from  the  increased  microearthquake 
coverage  during  the  1980’s  (Bungum  et  al. ,  1986),  including  the  most  re¬ 
cent  data  from  the  networks  and  stations  in  northern  Norway  (Kvamme  and 
Hansen,  1989). 

In  Figures  2.4.1-2.4.3,  there  are  some  patterns  of  seismicity  that  appear 
clearer  than  reported  before,  especially  when  compared  with  the  structural 
information.  In  addition  to  the  Mid-Atlantic  Ridge,  the  seismicity  follows 
broadly  the  continental  margin  from  Svalbard  to  the  northern  North  Sea. 
Some  concentrations  of  earthquakes  also  are  found  in  the  coastal  areas  of 
western  and  northern  Norway  as  well  as  in  the  Viking  Graben,  the  Central 
Graben  and  the  Oslo  Graben.  In  addition,  significant  earthquake  activity  is 
confined  to  a  large  area  in  the  eastern  Lofoten  Basin  west  of  the  Senja  Fracture 
Zone  and  possibly  also  to  a  small  area  near  the  intersection  between  the  V0ring 
Plateau  Escarpment  and  a  set  of  parallel  faults  tied  to  the  East  Jan  Mayen 
Fracture  Zone.  We  also  note  the  presence  of  relatively  well-defined  zones  in 
Sweden  (Slunga  et  al ,  1984)  and  in  Finland  (Ahjos  et  al ,  1984).  These  zones 
have  significantly  lower  recurrence  rates,  however. 

Because  of  the  particular  properties  of  the  data  bases  used  in  this  analysis, 
and  the  way  they  have  been  combined,  we  have  chosen  to  break  the  seismicity 
data  into  three  different  time  periods,  as  follows: 

-  Figure  2.4.I:  1880-1954,  consisting  mostly  of  the  felt  (historical)  earth¬ 
quakes,  with  a  few  additional  instrumental  network  locations. 

-  Figure  2.4-2:  1955-1979,  consisting  mostly  of  instrumentally  located 
earthquakes  from  a  variety  of  reporting  agencies. 

-  Figure  2.4-3:  1980-1989,  consisting  mostly  of  locations  based  on  new 
microearthquake  networks. 

The  earthquake  magnitudes  used  in  these  catalogues  are  mostly  Ms  for 
the  first  one,  mi,  for  the  second,  and  Mi  for  the  third.  The  Ms  scale  is  here 
tied  to  the  “felt  area”  Aiv  (for  MSK  intensity  IV)  and  Am  (for  intensity  III) 
as  developed  by  (Muir  Wood  and  Woo,  1987). 

The  following  relationship  has  been  developed  between  Ms  and  Mi  (Bungum 
and  Seines,  1988): 

Ms  =  0.85Ml  +  0.60  (1) 

This  magnitude  relationship,  developed  for  events  in  the  Ml  range  3  to  5, 
is  used  in  most  cases  when  earthquakes  from  different  catalogues  are  used 
together.  The  catalogues  also  contain  many  reported  body  wave  magnitudes 
me,  (Figure  2.2.2),  but  with  relationships  to  Mi  and  Ms  which  are  more 
uncertain. 
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Earthquake  location  precisions 

The  seismicity  information  presented  in  this  paper  covers  almost  200  years. 
In  the  beginning  of  this  time  period  Norway  was  sparsely  populated,  limiting 
the  possibilities  for  gathering  earthquake  information,  while  for  the  last  few 
years  the  country  has  been  very  well  covered  with  seismic  stations,  arrays  and 
networks.  The  earthquake  location  precisions  therefore  necessarily  also  vary 
considerably.  In  this  presentation  of  the  earthquake  catalogues  we  discuss 
three  different  time  periods,  and  the  main  reason  for  this  separation  is  that 
the  three  time  periods  also  represent  different  levels  of  precision  (as  well  as 
different  detectability  levels),  as  outlined  in  the  following: 

The  first  time  period  (1800-1954,  Figure  2.4.1)  consists  primarily  of  macro- 
seismically  determined  epicenters,  and  all  of  those  have  recently  (as  mentioned 
above)  been  reevaluated  in  a  consistent  way,  including  tuning  of  locations  us¬ 
ing  recent  instrumentally  located  (and  macroseismically  surveyed)  template 
earthquakes  (Muir  Wood  and  Woo,  1987).  Location  uncertainties  can  still  be 
quite  large,  however,  usually  in  the  range  20-60  km,  and  up  to  100  km  in 
particular  cases.  It  should  be  noted  here  that  large  (macroseismically  located) 
earthquakes  often  may  have  larger  uncertainties  than  smaller  ones.  For  the  few 
instrumentally  located  earthquakes  during  this  time  period,  the  uncertainties 
are  assessed  to  be  of  the  same  size  as  for  the  macroseismically  located  ones. 
This  catalogue  of  historical  earthquakes  has  been  found  to  be  complete  since 
about  1880  down  to  Ms  4.0  south  of  65°N  and  down  to  Ms  4.5  further  north 
(Muir  Wood  and  Woo,  1987). 

The  second  time  period  (1955-1979,  Figure  2.4.2)  contains  earthquake  re¬ 
ports  that  are  more  inhomogeneous,  consisting  of  locations  reported  by  a  vari¬ 
ety  of  seismological  agencies.  With  a  fairly  large  number  of  events  in  the  data 
base,  it  was  generally  not  possible  for  us  to  go  back  to  the  original  readings 
for  these  events.  The  sorting  and  screening  criteria  discussed  above,  however, 
were  designed  first  of  all  to  ensure  the  best  possible  quality  in  the  locations, 
which  partly  conflicts  with  the  criteria  that  would  provide  the  best  possible 
completeness.  The  resulting  location  uncertainties  are  estimated  to  be  in  the 
range  20-25  km  for  most  of  the  events,  but  also  less  than  10  km  in  some  cases. 

The  third  and  last  time  period  (1980-1989,  Figure  2.4.3)  contains  again 
data  that  are  more  homogeneous,  with  completeness  down  to  about  Mi  2.0 
for  the  areas  south  of  65°N  .  These  locations  are  based  on  recordings  from  the 
relatively  dense  microseismic  networks  installed  in  the  area  since  the  beginning 
of  the  1980s.  We  have  in  this  case  been  actively  involved  in  much  of  the 
earthquake  location  work  (e.g.,  Bungum  et  al ,  1986),  and  most  of  the  readings 
are  also  available  to  us.  This  makes  it  easier  to  assess  location  uncertainties, 
found  to  be  less  than  15  km  in  most  cases,  and  often  less  than  5  km. 

All  of  the  earthquake  locations  based  on  instrumental  data  are  obtained 
using  an  essentially  continental  crustal  model.  Biased  errors  are  still  not  found 
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to  be  significant,  however,  because  most  of  the  offshore  seismicity  is  close  to 
the  coast  (with  a  reasonable  azimuthal  coverage  of  stations),  and  because 
readings  from  stations  in  Great  Britain  (including  Scotland  and  the  Shetland 
Islands)  and  from  ocean-bottom  seismographs  routinely  are  used  in  the  loca¬ 
tions.  For  the  events  large  enough  to  provide  readings  at  teleseismic  distances, 
the  azimuthal  coverage  is  even  better. 

Regional  patterns  of  seismicity 

In  the  areas  south  of  63°N,  the  offshore  seismicity  is  more  or  less  confined 
to  the  shelf  edge  to  the  northwest  and  to  the  Viking  Graben  and  the  Central 
Graben  to  the  west  and  southwest.  In  the  Central  Graben,  the  earthquakes  are 
quite  small,  but  nevertheless  occur  along  a  convincing  lineation.  The  events 
in  the  Viking  Graben  are  larger  and  more  frequent,  including  a  widely  felt  Ms 
5.3  earthquake  in  1927.  Some  seismic  activity  is  also  found  in  the  Oslo  Rift. 
Presently  only  very  small  earthquakes  occur  there  (Bungum  and  Fyen,  1980), 
but  one  of  magnitude  Ms  5.4  occurred  on  October  23,  1904,  and  another  of 
Ms  4.4  on  March  6,  1953.  Another  seismically  active  area  is  west  of  Jutland 
in  the  Norwegian-Danish  Basin,  possibly  related  to  the  Fjerritslev-Tornquist 
Zone  or  the  Fennoscandian  Border  Zone.  The  areas  south  of  63°N  with  the 
highest  earthquake  activity,  however,  are  the  coastal  areas  between  59  and 
63°N,  including  an  offshore  extension  into  the  northern  North  Sea. 

The  area  between  63  and  70°N  includes  the  location  (near  66.5°N,  14.5°E) 
of  the  largest  known  (Ms  5.8  in  1819)  historical  earthquake  in  northwestern 
Europe  (Bungum  and  Seines,  1988).  It  is  seen  fro  the  seismicity  maps  that  this 
area  along  the  Nordland  coast  is  still  quite  active.  One  of  the  most  prominent 
features  is  the  zone  of  seismic  activity  that  runs  parallel  to  the  shelf  edge  from 
64°N,  6°E,  continuing  into  the  Lofoten  Islands. 

For  the  area  between  70  and  78°N,  the  most  obvious  feature  (the  spreading 
ridge  excluded)  is  the  seismicity  that  seems  to  be  limited  by  the  series  of  faults 
and  fracture  zones  that  parallels  the  shelf  edge  between  northern  Norway  and 
Svalbard.  The  most  prominent  of  these  is  the  Senja  Fracture  Zone,  along 
which  we  find  a  well-defined  seismic  activity  (see  also  Kvamme  and  Hansen, 
1989).  A  broad  zone  of  seismicity  in  the  eastern  Lofoten  Basin  (cf.  Husebye 
et  al ,  1975)  now  seems  to  be  even  more  clearly  defined  than  before.  Further 
north,  the  Svalbard  region  also  shows  significant  earthquake  activity  in  some 
areas  (Mitchell  et  al,  1989).  Further  south,  events  have  been  located  along  the 
Senja  Fracture  Zone  and  eastwards  to  the  Ringvasspy-Loppa  Fault  Complex. 
The  Senja  Fracture  Zone  is  known  to  be  quite  active,  separating  the  oceanic 
crust  to  the  west  from  the  continental  crust  to  the  east. 

For  a  more  in-depth  study  of  the  seismicity  and  seismotectonics  of  Norway 
and  surrounding  areas,  including  focal  mechanism  solutions  (stress  patterns) 
and  focal  depths,  we  refer  to  Bungum  et  al  (1991). 
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Conclusions 


The  new  historical  and  instrumental  earthquake  catalogues  presented  here 
have  been  organized  to  cover  the  time  periods  1880-1954, 1955-1979  and  1980- 
1989.  These  catalogues  have  given  us  a  much  improved  picture  of  the  distri¬ 
bution  of  earthquakes  both  in  time,  space  and  magnitude.  The  seismicity 
follows  on  a  regional  basis  various  geologic  structures  such  as  fault  zones  and 
complexes,  grabens,  fracture  zones,  shelf  and  margin  areas. 


H.  Bungum 
A.  Alsaker 
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Fig.  2.4.3.  Seismicity  of  Norway  and  surrounding  areas  for  tl  me  period 
1980-1989.  Magnitude  is  predominantly  Ml- 


2.5  Report  on  the  symposium  entitled  “Regional  Seismic  Ar¬ 
rays  and  Nuclear  Test  Ban  Verification”  held  in  Oslo, 
Norway,  14-17  February  1990 

During  14-17  February  1990  NORSAR  hosted  an  international  symposium 
entitled  “Regional  Seismic  Arrays  and  Nuclear  Test  Ban  Verification”.  The 
symposium  was  attended  by  76  scientists  and  representatives  from  21  coun¬ 
tries,  including  a  large  number  of  seismologists  participating  in  the  work  of  the 
Conference  on  Disarmament’s  Group  of  Scientific  Experts  (GSE)  in  Geneva. 

The  purpose  of  the  symposium  was  to  assess  the  state-of-the-art  of  research 
on  regional  seismic  arrays  and  associated  topics.  In  particular,  the  symposium 
focused  upon  the  advanced  regional  arrays  NORESS  and  ARCESS  in  Norway 
and  their  associated  data  processing  facilities,  in  the  light  of  the  potential  of 
such  arrays  to  provide  a  much  improved  monitoring  capability  for  a  future 
comprehensive  nuclear  test  ban  treaty.  During  a  three-day  scientific  sympo¬ 
sium,  a  number  of  presentations  were  given  on  topics  relevant  to  this  issue. 
A  special  session  was  devoted  to  summarizing  the  experience  and  discussing 
further  plans  for  the  on-going  international  GSE  experiment  (GSETT-2). 

In  this  paper,  we  give  a  brief  review  of  some  of  the  results  presented  during 
the  scientific  symposium.  The  majority  of  the  papers  presented  will  be  pub¬ 
lished  in  a  special  issue  of  the  Bulletin  of  the  Seismological  Society  of  America, 
scheduled  to  appear  in  December  1990.  Authors  and  titles  of  the  28  papers  of 
the  special  issue  are  given  in  the  Appendix. 

Development  of  regional  arrays 

Reviews  of  recent  developments  with  regard  to  regional  seismic  arrays  are 
presented  for  NORESS  and  ARCESS  in  Norway  [1],  GERESS  in  the  Federal 
Republic  of  Germany  [2]  and  FINESA  in  Finland  [3].  Paper  [1]  summarizes  the 
design  considerations  leading  to  the  establishment  of  the  first  regional  array, 
NORESS,  and  describes  how  the  success  of  this  new  array  concept  motivated 
the  deployment  of  additional  arrays  of  this  type.  The  paper  documents  the 
basic  signal  processing  techniques  used  in  real-time  data  analysis  for  regional 
arrays,  and  demonstrates  the  excellent  detection  performance  of  such  arrays 
at  regional  distances  (less  than  2000  km).  It  is  shown  that  NORESS  and 
ARCESS  are  capable  of  detecting  seismic  events  of  magnitude  2.5  with  90 
per  cent  probability,  if  these  events  occur  within  1000  km  distance,  whereas 
global  teleseismic  networks  have  much  higher  event  detection  thresholds.  The 
FINESA  array  is  also  documented  to  have  an  excellent  performance  [3],  and 
together,  these  three  arrays  are  capable  of  locating  weak  seismic  events  in 
Fennoscandia  very  accurately  (typically  to  within  10-20  km).  The  GERESS 
array  currently  under  development  shows  many  of  the  same  excellent  features 
[2],  and  will  contribute  further  to  an  excellent  regional  coverage  of  large  parts 
of  Northern  Europe. 
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Processing  of  data  from  a  network  of  regional  arrays 

Recent  technological  advances  have  allowed  very  sophisticated  processing  tech¬ 
niques  to  be  applied  in  detecting,  locating  and  identifying  seismic  events  using 
a  network  of  seismic  arrays  and  single  stations,  and  this  is  highlighted  by  the 
development  of  the  Intelligent  Monitoring  System  (IMS)  [4], [5].  Two  of  the 
goals  for  this  system  are  (1)  to  demonstrate  the  monitoring  performance  and 
capability  of  the  system  for  small  events  at  regional  distances  and  (2)  to  ex¬ 
plore  the  promise  of  an  expert-systems  approach  for  providing  improved  mon¬ 
itoring  performance  as  experience  accumulates.  The  first  operational  version, 
described  in  [4],  processes  data  from  NORESS  and  ARCESS,  whereas  later 
versions  will  be  expanded  to  networks  including  both  arrays  and  single  sta¬ 
tions.  The  IMS  is  ambitious  in  exploring  and  integrating  many  new  computer 
technologies,  and  the  validity  of  the  concept  is  documented  in  an  evaluation 
of  its  initial  operational  performance  [5]. 

Signal  analysis  methods 

A  number  of  presentations  addressed  methods  for  processing  seismic  signals 
recorded  by  arrays  as  well  as  three- component  stations.  It  was  demonstrated 
that  both  types  of  stations  can  provide  information  very  useful  in  phase  identi¬ 
fication,  azimuth  estimation  and  estimating  the  apparent  velocity  of  detected 
phases.  From  theoretical  considerations  as  well  as  from  experimental  com¬ 
parison  [12], (24], [14]  arrays  are  shown  to  be  superior  in  this  regard  at  low 
signal-to-noise  ratios,  although  the  precision  e.g.  of  azimuth  estimates  is  in¬ 
fluenced  by  a  number  of  factors,  including  phase  type,  frequency  of  the  signal 
and  systematic  bias  caused  by  earth  heterogeneities  [14], [17],  [24].  A  very 
promising  approach,  discussed  in  [16]  is  that  of  joint  analysis  of  3-component 
and  array  data. 

Signal  detection  methods  are  discussed  in  several  papers.  In  [ll],  a  system 
for  on-line  detection  and  signal  analysis  is  presented  as  applied  to  a  Soviet  3- 
component  station  in  Kazakhstan.  In  [13],  a  detection  technique  is  described 
using  NORESS  array  and  3-component  data.  A  statistical  approach,  using 
adaptive  techniques,  to  detection  processing  and  estimation  is  presented  in 
[7].  A  new  approach  to  obtain  precise  relative  location  estimates  of  seismic 
events,  using  high  frequency  recordings,  is  presented  in  [23]. 
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Source  identification 

Traditionally,  seismic  discrimination  research  has  focused  on  distinguishing 
between  earthquakes  and  underground  nuclear  explosions.  Under  a  Compre¬ 
hensive  Test  Ban  Treaty,  emphasis  will  be  on  detecting  and  identifying  weak 
seismic  events,  and  a  third  category,  large  chemical  explosions  for  industrial 
purposes  (e.g.  mining  work)  will  become  important  to  consider.  In  [10],  a 
very  promising  method  is  applied  to  NORESS  data  to  discriminate  between 
earthquakes  and  ripple-fired  quarry  blasts  (mining  events  consisting  of  several 
explosions  closely  grouped  in  space  and  time).  Using  spectral  characteristics  of 
the  signals,  an  “automatic”  discriminant  is  proposed  computing  the  likelihood 
that  ripple  firing  occurred  in  each  given  case. 

In  [8],  a  novel  approach  making  use  of  artificial  neural  networks  is  used  to 
develop  a  classification  procedure  between  earthquakes  and  mining  explosions. 
Also  in  this  approach,  the  spectral  characteristics  of  the  signals  form  the  basis 
for  the  discriminants.  The  neural  network  appears  to  improve  in  particular 
the  classification  of  outliers  in  the  population,  and  reduce  the  number  of  un¬ 
certain  events.  Application  of  neural  networks  in  improving  seismic  processing 
performance  is  also  addressed  in  [9]. 

Of  considerable  interest  for  source  identification  is  also  the  method  pro¬ 
posed  in  [15],  applying  transfer  functions  to  transform  e.g.  between  recordings 
of  presumed  single  explosions  and  ripple-fired  explosions,  and  also  between 
recordings  at  different  NORESS  sensors  for  a  given  event.  This  gives  promise 
to  improve  the  coherence  of  seismic  phases  recorded  at  an  array,  with  ensuing 
implications  for  improved  source  parameter  estimation.  In  [6],  a  case-based 
reasoning  approach  to  event  identification  is  discussed,  and  a  waveform  enve¬ 
lope  matching  technique  is  applied  to  a  set  of  Western  Norway  earthquakes 
and  explosions. 

Detection  thresholds  and  in-country  networks 

While  regional  arrays  were  originally  designed  to  enhance  the  capabilities  for 
detecting  and  characterizing  weak  seismic  events  at  regional  distances,  they 
have  also  been  found  very  effective  in  the  teleseismic  distance  range.  As  an 
example,  published  yields  of  Soviet  underground  nuclear  explosions  at  Semi- 
palatinsk  have  been  used  to  evaluate  the  NORESS  detection  threshold,  in 
terms  of  explosive  yield  for  events  at  this  test  site  [20].  The  threshold  for  de¬ 
tection  at  NORESS  is  estimated  to  be  as  low  as  0.1  kt,  assuming  full  coupling 
and  normal  noise  conditions.  It  is  pointed  out  that  NORESS  has  particularly 
favorable  conditions  for  detecting  small  events  from  this  test  site,  and  that  the 
seismic  identification  threshold  necessarily  will  be  higher  than  the  detection 
threshold. 

Data  from  new  Global  Seismic  Network  stations  in  the  Soviet  Union,  in¬ 
stalled  as  a  cooperative  project  between  American  and  Soviet  scientists,  have 
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been  applied  in  several  studies  to  address  problems  relevant  to  an  in-country 
monitoring  network.  Seismic  noise  levels  at  these  stations  are  analyzed  in 
[18],  and  found  to  be  higher  than  at  NORESS  in  the  band  1-20  Hz,  with  max¬ 
imum  difference  ranging  from  7  to  25  dB,  depending  on  the  station.  However, 
significant  noise  reduction  can  be  achieved  by  borehole  deployment. 

Using  data  from  stations  in  the  USSR,  the  frequency-dependent  attenua¬ 
tion  of  regional  seismic  phases  has  been  studied  in  [21].  Attenuation  charac¬ 
teristics  are  found  to  be  similar  to  those  observed  in  Scandinavia,  but  with  an 
absolute  Pn  amplitude  almost  a  factor  of  2  higher  in  eastern  Kazakhstan  for 
a  fixed  Lg  magnitude. 

Recordings  of  Semipalatinsk  nuclear  explosions  at  the  new  Global  Seismic 
Network  stations  in  the  Soviet  Union,  together  with  data  from  stations  in 
China  have  been  analyzed  in  [19]  and  it  is  shown  that  RMS  Lg  can  be  measured 
at  widely  separated  stations  with  a  remarkable  degree  of  consistency.  The 
standard  deviation  of  the  differences  between  pairs  of  stations  is  as  low  as 
0.03-0.04  in  logarithmic  units,  and  reliable  measurements  may  be  made  at 
magnitude  (m*,)  down  to  about  4.0  for  stations  situated  about  1500  km  away 
from  Semipalatinsk.  The  importance  of  this  observation  in  terms  of  supplying 
yield  estimates  for  nuclear  explosions  down  to  and  even  below  one  kiloton  is 
pointed  out. 

Earth  structure,  wave  propagation,  scattering 

Several  of  the  papers  were  devoted  to  studies  of  general  problems  in  seismology 
and  geophysics,  in  areas  relevant  to  the  seismic  monitoring  issue.  The  structure 
of  the  crust  and  upper  mantle  in  parts  of  northern  Eurasia  is  addressed  in 
papers  [22],  [25]  and  [27],  All  three  papers  are  specifically  making  use  of 
regional  array  data.  Seismic  wave  propagation  and  scattering  are  addressed 
in  a  number  of  papers,  e.g.  [13],  [24],  [26],  [27],  [28]. 

Conclusion 

The  Oslo  symposium  demonstrated  the  considerable  progress  in  the  field  of 
seismic  monitoring  during  recent  years.  It  particularly  highlighted  the  tech¬ 
nological  advances  in  seismic  instrumentation,  data  communication  and  com¬ 
puter  processing,  as  exemplified  by  the  development  of  advanced  regional  seis¬ 
mic  arrays  with  very  sophisticated  automatic  and  interactive  signal  process¬ 
ing  facilities.  The  presentations  at  the  scientific  symposium  show  that  these 
technological  advances  are  accompanied  by  considerable  scientific  progress,  al¬ 
though  much  work  remains  in  order  to  fully  exploit  the  potential  offered  by 
regional  arrays  in  a  seismic  monitoring  context. 


F.  Ringdal 
S.  Mykkeltveit 
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Appendix 


List  of  papers  given  during  the  1990  Oslo  Symposium  on  Regional  Seismic 
Arrays  and  Nuclear  Test  Ban  Verification,  to  appear  in  a  special  issue  of  The 
Bulletin  of  the  Seismological  Society  of  America 
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2.6  Summaries  of  Quarterly  Technical  Reports  Submitted 

During  the  current  fiscal  year  (FY90),  three  quarterly  technical  reports  were 
submitted  on  this  contract.  The  abstracts  of  these  papers  are  given  in  the 
following. 


2.6.1  The  stability  of  RMS  Lg  measurements  and  their  poten¬ 
tial  for  accurate  estimation  of  the  yields  of  Soviet  under¬ 
ground  nuclear  explosions  —  R.A.  Hansen1,  F.  Ringdal1  and 
P.G.  Richards2  NTNF/NORSAR,  Kjeller,  Norway,  and 
2  Lamont-Doherty  Geol.  Obs.,  Columbia  Univ.,  New  York, 
USA) 

Data  on  underground  nuclear  explosions  have  recently  become  available  from 
modern  digital  seismic  stations  installed  within  the  Soviet  Union  and  China. 
Observations  of  root  mean  square  (RMS)  Lg-wave  signals  for  Soviet  under¬ 
ground  nuclear  explosions  at  the  Shagan  River  test  site  in  East  Kazakhstan 
show  that  the  relative  amplitudes  of  the  RMS  signals,  at  stations  in  Norway, 
the  U.S.S.R.  and  China,  are  very  similar  for  different  explosions,  the  standard 
deviation  of  the  differences  being  only  about  0.03  in  logarithmic  units  (i.e., 
magnitude  units). 

This  is  consistent  with  earlier  observations  comparing  NORSAR  and  Grae- 
fenberg  array  data,  and  the  observed  scatter  is  significantly  lower  than  has 
been  reported  for  Lg  data  from  Nevada  Test  Site  explosions.  In  view  of  the 
excellent  correspondence  found  by  Nuttli  (1986)  and  Patton  (1988)  for  Lg 
versus  yield  at  Nevada,  this  indicates  that  RMS  Lg  has  a  potential  for  yield 
estimation  with  very  high  accuracy  at  Shagan  River. 

Our  study  has  shown  that:  (a)  selected  stations  in  the  U.S.S.R.  and  China, 
situated  at  regional  distances,  provide  a  much  improved  signal-to-noise  ratio  of 
the  Lg  phase  for  events  at  Shagan  River,  as  compared  to  NORSAR  array  data; 
(b)  the  scaling  of  RMS  Lg  amplitudes  between  different  sized  events  recorded 
at  the  same  single  station  site  appears  to  be  consistent  with  that  of  NORSAR, 
indicating  a  remarkable  degree  of  precision  in  single  station  measurements  of 
Lg  signal;  (c)  RMS  Lg  amplitude  measurements  for  the  best  of  these  stations 
may  be  made  at  1.5  to  2.0  magnitude  units  lower  than  at  NORSAR  or  Grae- 
fenberg,  allowing  a  much  lower  threshold  for  Lg  based  yield  determinations; 
and  (d)  the  P-wave  detection  capabilities  of  these  single  stations  do  not  match 
those  of  the  NORESS  and  ARCESS  arrays,  thus  teleseismic  signals  continue 
to  be  important  for  the  detection  of  small  nuclear  explosions. 

Our  conclusion  is  that  Lg  signals  appear  to  provide  an  excellent  basis  for 
supplying  estimates  of  the  yields  of  nuclear  explosions  even  down  to  below 
one  kiloton,  when  such  signals  are  recorded  at  high-quality  digital  in-country 
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seismic  stations,  and  when  calibrated  by  access  to  independent  (non-seisinic) 
yield  information  for  a  few  nuclear  explosions  at  the  test  sites  of  interest. 
In  the  context  of  monitoring  a  low  yield  threshold  test  ban  treaty,  it  will, 
in  addition,  be  important  to  take  into  consideration  various  environmental 
conditions  in  the  testing  area,  such  as  the  possible  presence  of  cavities,  and  to 
devise  appropriate  procedures  for  on-site  observations  in  this  regard. 


2.6,2  The  Mr  scale  in  Norway  —  A.  Alsaker,  L.B.  Kvamme,  R.A. 

Hansen,  A.  Dahleand  H.  Bungum  (NTNF/NORSAR,  Kjeller, 

Norway) 

A  new  low  magnitude  Mt  scale  has  been  developed  for  Norway,  based  on 
a  regression  analysis  of  synthesized  Wood-Anderson  records.  The  scale  is 
applicable  for  distances  up  to  more  than  1000  km,  and  the  data  used  comprise 
741  short-period  recordings  at  21  seismic  stations  from  195  earthquakes  in 
the  magnitude  range  1  to  5  occurring  in  and  around  Norway  over  the  last  20 
years.  Magnitude  corrections  for  distance  have  been  evaluated  in  terms  of  a 
geometrical  spreading  term  a  and  an  anelastic  attenuation  term  6,  and  the 
significant  regional  crustal  differences  in  the  area  under  investigation  made  it 
desirable  to  develop  these  for  several  subsets  of  the  data  base.  The  results  for 
a  are  generally  found  to  be  around  the  commonly  found  value  of  1.0  (using  the 
Lg  phase),  while  the  6-values  are  found  to  be  around  0.0008,  consistent  with 
the  weak,  intraplate  attenuation  expected  for  Norway.  Compared  to  interplate 
California,  this  difference  in  attenuation  represents  more  than  one  magnitude 
unit  at  a  distance  of  1000  km. 

New  Mj,  scales  are  commonly  tied  to  Richter’s  original  definition  at  the 
standard  reference  hypocentral  distance  of  100  km.  The  significantly  weaker 
Lg  wave  attenuation  in  Norway,  however,  requires  a  smaller  reference  distance. 
We  have  chosen  a  value  of  60  km,  based  on  an  overall  assessment  of  regional 
coverage,  focal  depths  and  quality  of  the  data.  The  resulting  M i  formula  for 
Norway  reads 

ML  =  logAwa  +  alog{R/ 60)  +  6  (.ft  -  60)  +  2.68  +  5  (1) 

where  Awa  is  synthesized  Wood-Anderson  amplitude  (in  mm),  R  is  hypocentral 
distance  (in  km),  and  5  is  a  station  correction  term  that  for  all  21  stations 
is  found  to  lie  within  the  range  ±  0.22.  When  using  the  entire  data  base 
the  spreading  term  a  equals  1.02,  and  the  anelastic  attenuation  term  6  equals 
0.00080.  When  only  strictly  continental  ray  paths  are  selected,  the  a-value 
decreases  to  0.91  while  the  6-value  increases  to  0.00087,  a  difference  which 
on  the  average  accounts  for  less  than  0.1  magnitude  units.  While  all  values 
used  in  the  regressions  have  been  derived  for  vertical  amplitudes,  a  separate 
analysis  has  shown  that  these  are  not  significantly  different  from  the  horizontal 
ones,  and  the  new  scale  is  therefore  applicable  to  both.  In  order  to  facilitate 
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the  practical  use  of  this  new  Ml  scale,  a  relation  has  also  been  established 
between  observed  seismogram  amplitudes  (corrected  for  instrument  response) 
and  the  synthesized  Wood- Anderson  amplitudes.  This  relation  reads  logAwa  — 
0.925  logAobs  -  2.32. 

The  new  Ml  scale  magnitudes  for  the  events  analyzed  are  in  fairly  good 
agreement  with  those  calculated  from  a  previously  used  relation  developed  by 
Bath  for  Sweden.  The  new  values  are,  however,  systematically  about  0.4  mag¬ 
nitude  units  lower,  which  is  mostly  due  to  the  combined  effect  of  a  reference 
distance  less  than  100  km  and  a  Wood- Anderson  magnification  of  2080  instead 
of  the  earlier  value  of  2800.  The  new  Ml  values  have  also  regressively  been 
related  to  a  data  set  of  M,  values,  yielding  the  relation  M ,  =  0.83 Ml  +  1.09. 


2.6.3  Application  of  regional  arrays  in  seismic  verification  research 

—  S.  Mykkeltveit1,  F.  Ringdal1,  T.  Kvaerna1  and  R.W.  Alewine2 
(i  NTNF/NORSAR,  Kjeller,  Norway,  and  2  DARPA/NMRO, 
Virginia,  USA) 

The  paper  gives  an  account  of  the  work  related  to  the  development  of  the 
NORESS  concept  of  a  regional  array.  The  array  design  considerations  and  ob¬ 
jectives  are  reviewed,  and  a  description  is  given  of  the  NORESS  and  ARCESS 
array  facilities  in  Norway  with  their  field  installations,  data  transmission  lines 
and  data  receiving  center  functions. 

The  automatic  signal  detection  processing  of  NORESS  data  applies  mul¬ 
tiple  narrow-band  frequency  filters  in  parallel  and  forms  array  beams  from 
selected  subgeometries.  The  detection  algorithm  is  based  on  computing  the 
STA/LTA  ratio  for  each  beam  individually,  and  a  detection  is  declared  when- 
i  r  this  ratio  exceeds  a  given  threshold.  It  is  explained  how  the  beam  de¬ 
ployment  and  the  individual  threshold  values  can  be  tuned  to  ensure  that  the 
interesting  phase  arrivals  are  not  missed,  but  at  the  same  time  avoiding  coda 
detections. 

For  each  detected  signal,  frequency- wavenumber  analysis  is  invoked  to  de¬ 
termine  arrival  azimuth  and  apparent  velocity.  Currently,  a  broad  band  esti¬ 
mator  is  used,  and  it  is  demonstrated  that  use  of  this  algorithm  increases  the 
stability  of  the  azimuth  and  apparent  velocity  estimates,  relative  to  narrow 
band  methods.  Local  and  regional  events  are  automatically  located  on  the 
basis  of  identification  and  association  of  P-  and  S-wave  arrivals.  The  uncer¬ 
tainty  in  the  arrival  azimuth  is  the  limiting  factor  in  accurately  determining 
single-array  event  locations,  and  it  is  shown  that  this  uncertainty  is  as  large 
as  10°-15°  for  Pn  phases  from  certain  regions. 

In  order  to  further  investigate  the  potential  of  the  NORESS  concept,  work 
was  initiated  towards  installing  a  network  of  regional  arrays  in  the  northern 
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Europe  area.  This  involved  the  deployment  of  the  ARCESS  array  in  northern 
Norway,  and  the  installation  of  the  FINESA  array  in  Finland  in  cooperation 
with  the  University  of  Helsinki.  Data  from  these  three  arrays  have  been  used 
jointly  in  a  location  estimation  scheme.  It  is  shown  that  for  events  in  the 
Fennoscandian  region  of  magnitude  typically  around  2.5  and  for  which  at 
least  one  phase  is  detected  by  each  array,  location  estimates  can  be  obtained 
automatically  that  deviate  from  published  network  locations  by  only  16  km 
on  the  average. 

In  the  future,  it  is  anticipated  that  additional  arrays  and  single  stations 
in  the  northern  Europe  area  will  contribute  real-time  data  to  NORSAR  for 
analysis  jointly  with  existing  arrays.  The  first  additional  data  to  become 
available  will  be  from  the  GERESS  array,  which  will  be  established  in  the 
Federal  Republic  of  Germany  in  1990.  Future  perspectives  also  include  the 
use  of  expert  system  technology  in  the  data  analysis,  and  the  IMS  system 
already  in  operation  represents  the  initial  attempt  in  this  regard.  A  summary 
is  given  of  problem  areas  where  further  work  is  needed  in  order  to  fully  exploit 
the  regional  array  concept. 
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